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1 . 0 SUMMARY 


1.1 Identification of Program 

This is the final report for the Analytical Design of an 
Advanced Radial Turbine study program under NASA Contract 
DEN3-106, to define an optimum ceramic radial inflow, single- 
stage turbine system, for an advanced automotive gas turbine 
engine. The turbine study program is based on an advanced gas 
turbine cycle, 100.67-kW (135-shp) engine for use in a 1588-kg 
(3500-lb) automobile, with a specific fuel consumption of 
0.227 kg/kW-hr (0.373 Ib/hp-hr) . 

Utilization of high turbine cycle temperatures of 1370®C 
(2500°F) and ceramic materials are prime features of this study. 

1 . 2 Objective 

The program objective is to investigate the effects of tip 
speed, non-radial rotor blading, inducer- to-exducer work split, 
and deswirl vanes on the system efficiency and mechanical reli- 
ability of a single-stage, ceramic, radial turbine designed to 
meet the performance requirements of an advanced automotive gas 
turbine cycle. 

1 . 3 Summary of Major Results 

Analysis results indicate that, based on projected 1983 
aerodynamic performance and ceramic material properties, radial 
blade rotor configurations with a tip speed of 701 m/sec 
(2300 fps) are feasible and satisfy the goals of 87.0-percent 
system efficiency and 0.9999 cumulative probability of success 
mechanically. 

Nonradial blade rotor configurations allow rotor inducer 
tip speed to be reduced for equivalent performance. However, 
the inherent higher blade stress associated with this type of 
configuration requires a drastic increase in blade thick- 
ness to restore reasonable cumulative probability of success 
(CPS) values for the ceramic material blade. The required 
increase in blade material thickness, in turn, results in higher 
disk loadings and reduced CPS values for the wheel. Therefore, 
from aerodynamic, mechanical, and off-design considerations, 
the radial-blade rotor appears to offer the lowest overall risk 
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and is the optimum configuration 
design includes the following: 

Shaft speed 

Inducer tip speed 

Rotor inlet blade angle 

Rotor blade number 

Rotor exit swirl 

Rotor exit critical velocity 

Rotor exducer-to-inducer tip 
radius ratio 


considered in this study. This 

\ 

95,000 rpm 
701 m/sec (2300 fps) 

0 degrees 
12 to 14 
-10 degrees 
ratio 0.449 
0.632 
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2 . 0 INTRODUCTION 


2. 1 Background 

The need for more efficient use of available fuel supplies 
has led to increased interest in gas turbine engines for auto- 
motive applications. Both government and industry have con- 
ducted studies to determine optimum cycle parameters and mechan- 
ical configurations. A high-temperature, regenerated, single- 
shaft engine has created considerable interest for automotive 
application due to economy, cost, packageability , and drive- 
ability. In a current NASA-DOE contract design study effort for 
an improved gas turbine engine, a single-stage radial turbine 
has emerged as the leading candidate compared to a single- or 
multi-stage axial turbine. 

The single-stage radial inflow turbine has attributes par- 
ticularly suited to this size engine because of ruggedness, 
lower sensitivity to tip clearance, and higher stage work levels 
when compared to axial turbines. 

For a single-shaft automotive application, a single-stage 
radial turbine must effectively cope with the basic requirement 
of high-stage work levels at maximum power and still maintain 
acceptable performance over the entire automotive duty cycle to 
achieve significant reductions in fuel consumption. However, by 
reducing fuel consumption on the order of 30 to 50 percent, the 
gas turbine will offer a distinct advantage in fuel economy over 
current spark-ignition engines and provide an automotive engine 
that is adaptable to a wide range of alternate fuels. 

2.2 Related Work 


For a single-shaft configuration, studies conducted to date 
show that work levels from 419 to 512 kJ/kg (180 to 220 Btu/lb) 
are required at maximum power for the radial turbine. Numerous 
experimental studies conducted at AiResearch have shown that 
maximum radial turbine performance correlates well with the cen- 
trifugal compressor slip factor developed by Stanitz (ref. 1) . 
That is, peak radial turbine performance occurs at a rotor inlet 
work coefficient (rotor inlet tangential velocity /rotor inlet 
tip speed) that is consistent with the slip factor criteria 
developed by Stanitz. On this basis, a relationship between 
stage work and inducer tip speed can be derived for maximum 
radial turbine efficiency. 

This relationship is presented in Figure 1 for zero exit 
swirl and for a range of inducer blade numbers. At the expected 
maximum power work levels, tip speeds between 671 and 762 m/sec 
(2200 and 2500 fps) would be required to achieve peak radial 
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kj/kg 0 233 279 326 372 419 465 512 558 

SPECIFIC STAGE WORK, AH-]- 

Figure 1. Optimum Tip Speed as a Function of 
Stage Work and Blade Number. 
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turbine performance. Cycle studies indicate that high turbine 
inlet temperatures [1370®C (2500®F)] also will be required to 
achieve the fuel reduction goals. Ceramic materials have 
received considerable attention for this application because 
they eliminate the need for turbine cooling and the critical 
materials required for metal turbines, as well as having the 
potential for low cost. However, at the start of the program, 
the mechanical properties of ceramic materials were in a period 
of such rapid improvement that an accurate estimate of maximum 
allowable tip speed of a radial turbine was not possible. 

If tip speeds from 549 to 610 m/sec (1800 to 2000 fps) were 
required to satisfy the material properties, severe tip speed 
limitations would exist relative to the optimum tip speed 
required to achieve maximum performance (Figure 1) . The impact 
of tip speed limitations on radial turbine performance is illus- 
trated in Figure 2. Results from a research turbine test 

program (designed with optimum tip speed) indicate that peak 
efficiency was achieved at the design point. However, increas- 
ing or decreasing tip speed from the design point results in a 
severe performance penalty. 

The radial turbine performance characteristics, shown in 
Figure 2, raise a number of questions: 

o What are the physical phenomena that cause the rapid 
decrease in performance at higher or lower than opti- 
mum rotor inducer tip speeds? 

o Once the physical phenomena are established, what 
design technique would tend to minimize the effect and 
reestablish maximum performance at nonoptimum tip 
speeds? 

o What technique would tend to optimize turbine perform- 
ance not only at maximum power but throughout the 
entire duty cycle range? 

These underlying questions formed the basis for arriving at 
a radial turbine design technique that would optimize radial 
turbine performance for the entire duty cycle of a single-shaft 
automotive gas turbine engine. 

First, assume that the effects of lower or higher than 
optimum tip speed are attributed to rotor inlet incidence loss, 
and then, that the incidence penalty is directly proportional to 
the difference in kinetic energy between the optimum inlet work 
coefficient and that actually implied by the rotor inlet vector 
diagram. Figure 2 shows that, with these assumptions, the 
predicted and measured decrement in turbine efficiency are in 
good agreement (considering that the effect of stator loss and 
rotor reaction are also present under these conditions) . Given 
this correlation between nonoptimura rotor tip speed and rotor 
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Figure 2. Effect of Rotor Tip Speed on Attainable Turbine Efficiency 



inlet incidence loss, what design technique would minimize these 
effects? Two approaches appear to be feasible: First, the 
rotor inducer-to-exducer work split could be adjusted to the 
point that optimum rotor inlet conditions would be reestablished 
(i.e., high negative rotor exit swirl). Second, the rotor inlet 
blade angle could be adjusted to reestablish optimum rotor inlet 
conditions with a 90-degree radial inlet flow configuration. 

Adjusting the inducer-to-exducer work split should be bene- 
ficial for all types of rotor blading as long as a limitation on 
tip speed or blade angle exists. Unfortunately, numerous exper- 
imental studies indicate that the turbine downstream ducting is 
relatively sensitive to rotor exit swirl, and that the magnitude 
of duct loss is highly dependent on the type of geometry uti- 
lized. For example. Figure 3 illustrates the loss characteris- 
tics as a function of inlet swirl for an interstage duct and 
exhaust diffuser. Although Dovzhik's data (ref. 2) are for uni- 
form inlet conditions, rig test results with the actual rotor 
inlet conditions show the same basic characteristics. Unpub- 
lished AiResearch results for radial diffusers indicate that 
this type of exhaust duct could be the least sensitive to inlet 
swirl; however, additional analyses and experimental studies are 
required to further define the design techniques. With the 
assumption that the turbine exit ducting loss is a function of 
swirl, the duct loss characteristics must be included in the 
optimization process for tip speed-limited radial rotor designs. 
Even if optimum tip speed and zero exit swirl are achieved for a 
given design point, the wide range of operating conditions 
experienced by the single-shaft turbine would still result in 
relatively high rotor exit swirl during the normal duty cycle. 

A turbine system optimization technique based on combining 
the performance characteristics of the turbine stage and exhaust 
duct is illustrated in Figure 4. At the corrected work level 
shown, a tip speed of 781 m/sec (2562.0 fps) would be required 
to achieve optimum rotor inlet conditions and zero exit swirl. 
However, with an allowable tip speed of 610 m/sec (2000 fps), 
zero exit swirl would occur at a rotor inlet work coefficient of 
1.37. Although the exhaust duct loss would be minimized under 
these conditions. Figure 4 shows that stage efficiency rapidly 
decreases as this level of inlet work coefficient is approached. 
By combining the performance characteristics of the turbine 
stage and exhaust duct, the rotor inlet incidence effects are 
minimized, and system efficiency is improved by 4.0 points. 

Figure 4 was based on conventional radial element rotor 
blades. Further improvement in system efficiency is expected 
with three-dimensional (nonradial) blades. AiResearch has 
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Figure 3. Illustration of Swirl Effect on 
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Figure 4. Turbine System Optimization Technique 
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investigated the aerodynamic and mechanical feasibility of non- 
radial rotors, and the results indicate the following: 

o With careful control of the rotor blade turning dis- 
tributed through the rotor, rig test results show that 
the increase in performance is in good agreement with 
that predicted from a decrease in inducer incidence 
loss 

o Minimized rotor blade bending stress requires exten- 
sive rotor blade thickness optimization and blade 
restacking. At the present time, 20 degrees of rotor 
inlet blade angle appears to be the upper limit in 
terms of rotor mechanical design 

The recent study conducted on deswirl vanes by Mitchell and 
Soileau (ref. 3) indicates that duct losses could be signifi- 
cantly reduced for high exit swirl and Mach number turbine 
designs. Therefore, deswirl vanes will be included in the 
system optimization. 

2.3 Scope/Purpose 

The program objective was to investigate the performance 
potential of a high-work, ceramic, single-stage, radial turbine 
with radial and nonradial rotor blading, with and without 
deswirl vanes, based on a system optimization approach applic- 
able to a single-shaft automotive engine. This was accomplished 
with a parametric study by examination of the aerodynamic and 
mechanical effects of the following parameters: 

o A range of shaft speeds to evaluate rotor size effects 

o A range of inducer tip speeds to determine the effects 
on efficiency and mechanical integrity 

o A range of rotor inlet blade angles to determine the 
effects on inducer blade loading and the impact on 
rotor inducer blade stresses 

o A range of rotor inducer-to-exducer work splits to 

determine the effects of rotor exit Mach number and 
exit swirl on exhaust diffuser losses 

o A range of rotor inducer tip-to-exducer tip radius 

ratios for each value of shaft and tip speed 

o An examination of blade number effects 

To evaluate turbine performance effects only, all cycle 
parameters, with the exception of radial turbine efficiency and 
exhaust duct pressure loss, were held constant for the analysis. 


10 



3.0 DISCUSSION OF RESULTS 

3.1 Single-Shaft Engine Cycle Selection 

The maximum power point cycle selected for the advanced 
radial turbine is representative of operating conditions pro- 
jected for a gas turbine automotive engine produced in the 
mid-1980s for a 1588-kg (3500-lb) automobile. 

The engine specific fuel consumption (SFC) is 0.1965 kg/kW- 
hr (0.323 Ib/hp-hr) , and the net power is 101 kW (135 hp) . The 
turbine system design requirements and performance goals at this 
condition are presented in Table I. Station nomenclature used 
in Table I is defined in Figure 5. 

3 . 2 Configuration Definition 

The selected system configuration consists of a turbine 
inlet scroll, radial stage and downstream axial exhaust diffuser 
(with and without deswirl vanes). Alternate turbine inlet and 
exhaust configurations are possible; indeed, further automobile 
engine analysis may indicate these are required. However, the 
primary objective of this program is to evaluate the performance 
potential of a ceramic radial turbine through use of a system 
optimization procedur,e that could be used to evaluate modifica- 
tions that may be required by further automobile engine 
analysis . 

The effect of selecting either a scroll or an annular inlet 
is expected to be small for the present program since, in either 
case, the objective is to maintain a relatively low Mach number. 
Therefore, pressure loss with either type of inlet would be 
expected on the order of 1/2 to 1 percent. The type of exhaust 
diffuser configuration selected is more significant due to the 
difference in the levels of losses expected from an axial or a 
radial diffuser. Since the type of exhaust diffuser selected 
for the analysis will limit the generality of the results, both 
at the selected maximum power design point and off-design condi- 
tions, the selection will be based on available diffuser data 
with actual radial rotor exit conditions imposed. The selected 
exhaust duct configuration shown schematically in Figure 5 is a 
conical diffuser and is the type normally used for small auxili- 
ary power units (APUs) . Detailed test results are available for 
this type of diffuser. 

3.3 Optimization Study Parameter Range 

The design parameters evaluated in the optimization study, 
and the associated range, are presented in Table II. The lower 
limit of turbine shaft speed (80,000 rpm) corresponds to the 
lowest specific speed at which peak efficiency could be main- 
tained, based on available data. The upper limit of shaft speed 
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TABLE I. CERAMIC RADIAL INFLOW TURBINE REFERENCE 
CYCLE REQUIREMENTS 


Parameters 

Units 

Requirements 

Turbine Operating Life 



hr 

3500 

Turbine Design 

Reliability 



0.9999 

Turbine Scroll 

Inlet Total 


°C ("F) 

1370 (2500) 

Temperature 






Cycle Pressure 

Ratio 




5:1 

Scroll 

Inlet Total Pressure 


kN/m^ (psia) 

478 (69.331) 

Rotor Exit Total Pressure 


kN/m^ (psia) 

116 (16.8) 

Mass Flow, w 




kg/sec (Ib/sec) 

373 (0.823) 

Specific Stage 

Work, Ah^ 



kJ/kg (Btu/lb) 

509 (218.73) 



1-^ 




Stage Total-to- 

-Total Pressure 


4.087 

Ratio, P_ 







T-T 

1-4 





Stage Inlet Corrected Flow, 


kg/sec (Ib/sec) 

0.191 (0.421) 

w/0/6 2 _ 





Stage Corrected Work, AH^ 

’/e 

1-4 

kJ/kg (Btu/lb) 

89.162 (38.333) 

Stage Total-to- 

■Total Efficiency, 


0.897* 

^T“T 

1-4 






Scroll 

Total Pressure Loss, 



0.01* 

AP/P 

0-1 






Exhaust Diffuser Total Pressure 


0.04* 

Loss, AP/Pj 4 _g (with or with- 
out deswiri vanes 



System 

Total-to-Total Efficiency 


0.870** 

TlT_ip 

0-6 






*Varies in parametric study 



|**Goal 

from parametric study 
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Figure 5 


Definition of Radial Turbine System for Parametric Study 


TABLE II. PARAMETER RANGE FOR OPTIMIZATION STUDY 


Parameters 

Units 

Range 

Rotational Speed, N 

rpm 

80,000 to 110,000 

Inducer Tip Speed, U. 

m/sec 

(fps) 

792 (2600) 

(optimum with 10 blades) 




753 (2470) 

(optimum with 20 blades) 




549 (1800) 
(lower limit) 

Inducer Blade Angle/ |3^ 

deg 

0 to 20 

Rotor Exit Swirl, a 

deg 

0 to -50* 

Rotor Exit Absolute 
Critical Velocity 
Ratio, V/A^^ 



0.30 to 0.60* 

Rotor Blade Number, N^ 



10 to 16 

Rotor Exducer-to-Inducer 
Radius Ratio, R^ /R„ 

4 -^3 



Minimum based on 50% 
rotor exit hub blockage, 
maximum of 0.80 

Axial Shroud Clearance 
(Ca) 

cm 

(in. ) 

0.0254 (0.010) 
(constant) 

Radial Shroud Clearance 
(Cr) 



0.0254 (0.010) 
(constant) 

Rotor Backface Clearance 
(Cb) 

1 


0.0508 (0.020) 
(constant) 

*Estimated; absolute range will result 
coefficient and annular area imposed. 

from the rotor exit work 
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(110,000 rpm) is based on maintaining acceptable single-stage 
centrifugal compressor design characteristics. Higher speeds 
are possible; however, this would involve a compressor perfor- 
mance trade-off (due to the increasing inducer tip Mach number, 
which is considered beyond the scope of the current program) . 
The range of selected shaft speeds is indicated on the specific 
speed curve in Figure 6. 

3.4 Aerodynamic Performance Correlations for Parametric Study 

The initial parametric study was based on a one-dimensional 
aerodynamic performance evaluation model for the turbine scroll, 
stage and axial exhaust diffuser system. Correlations for the 
scroll and exhaust diffuser were incorporated into an existing 
radial turbine design program to arrive at a system performance 
for the range of parameters presented in Table II and based on 
the selected cycle conditions shown in Table I. The correla- 
tions and methods utilized for the parametric study are; 


o 


Scroll Loss ; The pressure loss in a turbine scroll is 
primarily a result of turning the flow and wall fric- 
tion. Rogo (ref. 4) presented a method for predicting 
the turning loss for a turbine scroll based on a mass 
averaged dynamic head within the scroll. Eckert 
(ref. 5) developed a technique for predicting the 
friction loss for a compressor scroll. A technique 
for predicting total turbine scroll loss has been 
derived from the work of Eckert and Rogo in terms of 
a scroll inlet Mach number, scroll inlet diameter 
(a) at the tongue, scroll inside radius (r©) / and fric- 
tion resistance coefficient (A), as follows; 


AP 

’total 


where; 


c = In 


INLET ro'^/ 

(^) 


0.266 + 


2ttX 


(-0 - I ) 


Xt, = 


Xr. = 


[0.8686 ln(j) + 1.74] 
0.316 


(Rough) 


(Smooth) 


and 


Re = 
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<Ti 



RPM (ft3/4/5egl/2)-^ 30 40 50 60 70 80 90 100 110 

0.233 0.310 0.388 0.465 0.543 0.620 0.698 0.775 0.853 

SPECIFIC SPEED 


Figure 6. Specific Speed Correlation for Radial Turbines. 




The loss predicted from the preceding equation would 
be the minimum loss that would be expected for a 
scroll of circular cross section and design, using the 
continuity and radial equilibrium scroll design equa- 
tions. 

o Maximum attainable (base) efficiency is established 
from specific speed correlations based on NASA 
(ref. 6) and AiResearch data. The NASA correlation is 
presented in Figure 6. 

o Base efficiency is then corrected for Reynolds number 
effects based on results from NASA (ref. 7) and 
AiResearch data. The correlation is presented in Fig- 
ure 7. 

o Rotor inlet incidence loss is calculated for a speci- 
fied rotor inlet work coefficient from the following 
relationships: 

In general, 


X 


3,i 


S.F. 


U. + tan 
^3 ^3 


B 



For radial blades. 


^3,i 


= S.F. 


= 1 - — 
Nr 


(Stanitz Slip Factor) 


^3, act 




n 


i 


1 


1 + 




X 



2 
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Figure 7. Reynolds Number Correlation for Radial Turbines 




where: 


U = Rotor inducer tip speed 

Vp = Rotor inlet absolute radial velocity 

Vu^ = Rotor inlet absolute tangential velocity 

R- = Rotor inlet tip radius 
^3 

R^j^ = Rotor exit mean radius 

= Stage total specific work 

X 3 i = Rotor inlet work coefficient for zero 
' incidence loss 


o 


X 3 = Rotor inlet work coefficient from vector 
' diagram for a specific tip speed and 

inducer-to-exducer work split 


^ STG 

^4,M 


'f) 


i 




base 




B 


V 


R 


S.F. 


N 


B 


Stage work coefficient 

Rotor exit mean work coefficient 

Stage total efficiency with incidence 
loss 

Stage total efficiency before incidence 
loss 

Rotor inlet blade angle 

Rotor inlet radial velocity 

Slip factor for radial blades, 1 - (2/Ng) 

Number of rotor blades 


The effect of rotor shroud clearance is based on cor- 
relations derived from references 8 and 9 and on 
AiResearch test results. The performance penalties 
are a function of both axial (Ca) and radial (Cr) 
clearance. The experimental data from references 8 
and 9 are presented in Figures 8 a and 8 b. Additional 
rotor clearance effects are present with rotor scal- 
lops. The performance effects as a function of back- 
face clearance recently have been evaluated for the 
Model GTP305-2 turbine (ref. 10) and are presented in 
Figure 8 . The size and scallop depth of this turbine 
are similar to the turbine considered in this study. 


19 




RADIAL CLEARANCE. 

PERCENT OF EXDUCER PASSAGE HEIGHT 


(A) Summary of Radial and Axial Clearance Effects on Total Stage Efficiency (Ref. 8) 
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(B) Summary of Radial and Axial Clearance Effects on Total Stage Efficiency 
at Design Equivalent Values of Speed and Pressure Ratio (Ref. 9) 
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BACKFACE CLEARANCE, 
PERCENT OF INDUCER PASSAGE HEIGHT 


(C) Effect of Rotor Backface Clearance on Stage Total Efficiency 
at Design Conditions (Ref. 10) 


Figure 8. Summary of Radial Turbine Rotor Clearance Effects 
on Total Stage Efficiency. 
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In addition, the results of a 1978 AiResearch- 
sponsored radial turbine research program investi- 
gated the effects of scallop depth on radial turbine 
performance. These data will be used to supplement 
the correlation derived from the Model GTP305-2 tur- 
bine. 

o The performance effects of blade number are based on 
the Models GTCP305-1, GTCP36-4, and published Pratt 
and Whitney data (ref. 11). The resultant change in 
turbine efficiency represents the effect of increased 
blade loading as the blade number is reduced for a 
given work requirement. 

o Disk friction losses along the rotor backface are 
estimated based on the results of the investigation 
conducted in reference 12. This work is fully 
described in reference 13. Note, however, that disk 
friction is based on a full rotor backface disk. The 
effect of rotor scallops on rotor disk friction is 
currently not available. 

o The exhaust diffuser loss correlation was based on the 
Model GTC36-200 configuration shown in Figure 9a. The 
annular inlet to the diffuser is followed by a rotor 
exit hub dump, a constant area section, and a conical 
diffuser section. The constant area section down- 
stream from the rotor exit dump allows the flow to mix 
and stabilize, prior to the rapid diffusion near the 
diffuser exit. The net result is a linear static 
pressure rise along the diffuser shroud from rotor 
exit to diffuser exit. 

The non-dimens ionalized diffuser geometry and perfor- 
mance potential of this system is presented in 
Table III. The loss characteristics of the Model 
GTC36-200 exhaust diffuser as a function of inlet 
swirl are presented in Figure 10. Although a number 
of data points at high corrected speed or low-pressure 
ratios result in relatively low loss levels, the pre- 
dominant characteristic is at a high loss level, as 
shown. 

The lower loss levels would be equivalent to lower 

diffuser inlet blockage levels, which is consistent 

with diffuser test data. However, for the high work 

and pressure ratio levels of this program, the higher 

loss level is_considered realistic. The minimum loss 

coefficient (w . ) derived from Figure 10 is 0.30. 
min 
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COMBUSTION 

DISCHARGE 


COMBUSTION 

DISCHARGE 


1.865 

A4 


L = 
AR 


5.35 



— <L 



(B) 


RADIAL TURBINE SYSTEM PERFORMANCE FOR PARAMETRIC STUDY 
IS DEFINED FROM SCROLL INLET {COMBUSTION DISCHARGE) TO 
EXHAUST DIFFUSER EXIT PLANES (A) WITHOUT DESWIRL VANES' 

(B) WITH DESWIRL VANES. 

Figure 9. Exhaust Diffuser Configuration used in Parametric Study 






TABLE III. PERFOEMANCE AND NONDIMENS lONAL GEOMETRY 
FOR SELECTED DIFFUSER CONDIFURATION 


Parameters 


Value 


Exit Area/Inlet Area 


1.865 


L/AR 


Diffuser Length 
Inlet Height 


5.35 


Static Pressure Recovery at 
Design Conditions, R^ 


0.468* 


Maximum Static Pressure Recovery 


0.60** 


Diffuser Inlet Critical Velocity 
Ratio at Design-Point, V/A^^ 


0.463 


Average Diffuser Inlet Swirl at 
Design-Point (a) , deg 


-3.0 


*100-percent corrected Speed, Pr,j,_,P ^ ^ = 4.07 
'*110-percent corrected speed, Priji_ijiL ^ “ 3.85 










DUCT INLET SWIRL 


DEGREES 


Figure 10. Diffuser Duct Loss Data 






Mpj - AVERAGE INLET MACH NUMBER 


Figure 11. Deswirl Vane Loss Characteristics (ref. 3) 
Using Arbitrary Air-Foil Shapes. 






o 


when a rotor exit deswirl vane is used, the diffuser 
non-dimens ionalized parameters shown in Table III 
will be applied at the deswirl vane exit. The losses 
required for the deswirl vane to reestablish axial 
flow are based on the performance potential achievable 
with arbitrary airfoil shapes indicated in Figure 11. 
The test data shown are based on the work of Pratt and 
Whitney (ref. 3). The deswirl vane geometry (solid- 
ity, aspect ratio) and diffusion limits will be 
derived from current axial compressor vane tech- 
nology. A detailed description of the model used in 
the parametric study is presented in Appendix B. 

The above correlations represent the majority of the design 
data available to the radial turbine designer. The current 
state-of-the-art for radial turbine design is a result of two 
major factors. First, the number of design variables associated 
with the radial turbine are increased significantly compared to 
the axial turbine. For example, rotor clearance effects are a 
function of rotor inducer blade width, exducer blade height, 
backface seal clearance, scallop depth and scallop saddle con- 
figuration. As a result, extensive (and expensive) rig test 
programs are required to arrive at meaningful correlations, 
which tend to restrict the dissemination of radial turbine 
data. Secondly, the flow field in the radial turbine is highly 
three-dimensional, radial stators are usually in the very low 
aspect ratio region, and the relative vorticity complicates 
rotor inlet and exit flow predictions. 

Nevertheless, experience gained from axial turbine designs 
has shown that the effects of rotor exit blockage and stage 
reaction should not be ignored during the turbine design pro- 
cess. The effects of rotor blockage on overall radial turbine 
performance do not appear to be available in the open litera- 
ture. The rotor exit hub blockage of AiResearch turbines used for 
specific speed correlation is on the order of 30 percent, and 
under this condition, excellent performance levels are still 
achieved. However, for the automotive engine, turbine rotor 
weight and inertia will have a significant influence on engine 
acceleration and subsequent driveability; for this reason, a 
minimized rotor exit hub radius (lower rotor exducer hub-to- 
inducer radius ratio) is desirable. The resultant increase in 
rotor exit blade height would reduce the sensitivity to rotor 
radial clearance. 

Under these conditions, an order of magnitude effect for 
rotor trailing edge blockage is necessary to arrive at a 
rational tradeoff. Analytical analysis in terms of evaluating 
boundary layer characteristics along the rotor meridional 
streamlines is considered ineffective at the present time, due 
to highly three-dimensional flow characteristics in the rotor. 
The problem is further compounded by the local inducer flow sep- 
aration associated with the highly loaded designs. Even if this 
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approach were pursued, the time required to analyze a large num- 
ber of solutions on this basis would be prohibitive. Therefore, 
a simplified approach is taken, whereby the relative effects of 
rotor exit blockage are calculated, based on the compressible 
mixing loss due to the rotor trailing edge thickness and local 
relative velocity. The total pressure loss calculated in this 
manner is compared with the loss calculated for an existing ref- 
erence design, and the difference in calculated loss is used to 
penalize the higher blockage designs. The model will be based 
on modifying Stewart's (ref. 14) compressible mixing loss equa- 
tions to account for trailing edge thickness alone. 

The basis for this model is a result of numerous observa- 
tions of test result that indicate the boundary layer thickness 
in the rotor exit hub region is small. For example. Figure 12 
shows the local efficiency of the selected reference turbine 
derived from rotor exit survey data. As indicated, the effi- 
ciency is relatively constant at 95.5 percent up to approxi- 
mately 40 percent of the blade height and then decreases rapidly 
above this point. The implication from rotor powder traces and 
predicted hub velocity distributions is that high secondary flow 
generated at the rotor inducer region propagatfes to the rotor 
exit shroud region. The flow along the hub line then expe- 
riences high acceleration downstream from the rotor inducer 
region at a relatively low loading level. The net result is 
that the losses in the hub region are extremely low and account 
for the high efficiencies measured. 

On this basis, then, a calculated mixing loss due to the 
trailing edge thickness alone would be expected to provide a 
realistic estimate of the effects of increased blockage in this 
region. However, even with this simplified approach, the calcu- 
lation procedure becomes fairly complex in order to arrive at an 
estimated performance decrement for rotor exit blockage. For 
example, the rotor exit flow path dimensions, specified exit 
work coefficient, and continuity must be solved to satisfy the 
radial equilibrium equation inside the rotor trailing edge. 
However, a radial thickness distribution is required that satis- 
fies both the maximum tip thickness and maximum allowable blade 
stress with a radial distribution of blade angle (usually radial 
blades in this region) . Then, the calculated radial relative 
velocity in conjunction with the established radial thickness 
distribution can be utilized to calculate a relative total pres- 
sure loss from hub to tip based on blade blockage alone. 

Since this type of analysis could not be achieved with the 
one-dimensional performance model available for the parametric 
study and the effect of rotor blockage should not be ignored, 
the following procedure was adopted: 

o First, the one-dimensional model was used to establish 
peak performance characteristics for a range of tip 
speeds and rotational speeds. 
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ROTOR EXIT PASSAGE HEIGHT - PERCENT 



Vj.j, STAGE LOCAL TOTAL EFFICIENCY 


Figure 12. Representative Rotor Exit Efficiency Characteristics 
for Radial Turbines. 



o viable configurations were examined mechanically with 
two-dimensional finite element methods that estab- 
lished rotor blade thickness satisfying the maximum 
allowable blade stress. A separate computer program 
being developed was utilized to solve the radial equi- 
librium equation with the specified blade thickness, 
and an estimated rotor exit blockage performance 
decrement was calculated. 

The effect of stage reaction is expected to increase 
directly with increasing work coefficient. However, isolating 
this effect with highly loaded radial turbines is complicated by 
rotor inducer incidence and stator loss, which will change as 
inducer-to-exducer work split is adjusted to increase stage 
reaction. If a systematic procedure is used to extract known 
losses at tested off-design points, a relative rotor reaction 
effect should evolve. This approach has been applied to the 
Model GTP36-51 and 1976 research radial turbine test data to 
derive a preliminary correlation of stage reaction on radial 
turbine performance. Figure 13 shows results as desirable mean- 
line reactions of 50 to 60 percent. The stage reaction shown in 
Figure 13 is defined in the usual manner by the ratio of static 
enthalpy across the rotor, relative to total stage work: 

R = 

STG 

In terms of a stage work coefficient and relative velocity, the 
stage reaction can be expressed as: 

The effect of stage reaction was incorporated in the performance 
model and the results of the parametric study will be presented 
with and without these effects. 

An additional effect that cannot be accounted for is rotor 
shroud curvature. The shroud curvature, as a function of rotor 
size and exducer tip-to-exducer tip radius ratio, increased with 
increasing radius ratio and resulted in higher velocity peaks 
and subsequent diffusion on the suction surface shroud stream- 
line, as the radius ratio approached 1.0. Therefore, if the 
parametric study indicates that high radius ratio designs are 
more favorable, this effect would not become evident until the 
detail rotor flow solution was in progress. To minimize this 
occurrence, the radius ratio was limited to 0.80. 


STG 


2X 


STG 
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3. 5 Projected 1983-1985 Performance Improvements 


The previously discussed performance model is based on cur- 
rent 1979 technology. However, as indicated in Table I, the 
high-stage work and performance level required by the turbine to 
achieve a 50-percent reduction in fuel consumption is aggressive 
and will require significant improvements in turbine efficiency. 
With the smaller size associated with the automotive applica- 
tion, both stator aspect ratio and rotor tip clearance penalties 
will tend to restrict the attainable efficiency level, even if 
optimum tip speed is achieved. 

Historically, radial nozzle designs have been based on con- 
stant section profiles with parallel end-walls. These condi- 
tions, coupled with the favorable nozzle inlet-to-exit area 
ratio, result in a theoretically high-acceleration two- 
dimensional nozzle. However, in practice, the nozzle aspect 
ratio is extremely low compared with an axial turbine, and the 
rotor shroud curvature imposes cross-passage static pressure 
gradients at the stator exit and transforms the highly acceler- 
ating two-dimensional flow to a complex three-dimensional flow 
f ield. 

Indications that the stator is not the low- loss two- 
dimensional device originally expected are implied from flow 
visualization techniques (powder traces) and from an inability 
to match radial turbine off-design test data utilizing tech- 
niques similar to reference 15, unless stator loss coefficients 
are increased significantly. Powder traces from a 0.30 aspect 
ratio radial nozzle showed that secondary flow accumulations 
covered approximately 90 percent of the passage in the trailing 
edge suction surface region. Even with constant section pro- 
files and parallel end-walls, secondary flows would be expected 
due to the end-wall pressure-to-suction static pressure grad- 
ients that result from the high turning. In addition, studies 
conducted by Langston (ref. 16) and Gaugler (ref. 17) show that 
blade-to-blade secondary flows are further compounded by the so- 
called horseshoe vortex generated at the nozzle leading edge. 

The concept of stator end-wall contouring and non-free 
vortex work distributions has been widely utilized in axial tur- 
bine designs over the past few years with excellent success. 
Since the effect of end-wall contouring is to reduce the blade- 
to-blade static pressure gradient (loading) at the nozzle end- 
walls, this concept will apply to radial nozzles also. Indeed, 
the benefits may be more pronounced with radial nozzles since 
the aspect ratios are significantly lower compared with axial 
nozzles. Therefore, in 1978 the potential benefits of end-wall 
contouring were investigated for radial stators with a compress- 
ible 3-D viscous flow analysis computer program developed by 
AiResearch under USAF Contract (F33615-76-C-2110) . The study 
was conducted as part of the TARADCOM advanced radial turbine 
program (ref. 18). Five different end-wall configurations were 
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examined (including a parallel wall baseline case) and compared 
on a predicted total loss basis. The analytical results showed 
that total pressure loss is reduced by 50 percent with end-wall 
contouring. Additional analytical and experimental effort is 
required to identify the true potential of the concept and also 
evaluate the performance potential of combining both the end- 
wall contouring and three-dimensional vane profile (non-free 
vortex type) concepts. However, on the basis of preliminary 
analysis, a 1.0-point improvement in turbine performance is con- 
sidered feasible. 

Although the sensitivity to rotor shroud clearance is 
reduced compared with axial turbines, scaling the allowable run- 
ning clearance with the reduced size of the automotive turbine 
may not be possible. In addition, the rotor scallop introduces 
a third clearance effect that is a function of both backface 
seal clearance scallop depth and scallop saddle configuration. 

Under these conditions, the estimated clearances shown in 
Table II will result in relatively high performance decrements. 
If higher engine running clearances are required due to shaft 
dynamics or differential thermo growth, then attainable radial 
turbine performance will be appreciably effected. Therefore, 
the concept of clearance treatment to minimize clearance effects 
will be considered for this application. 

Dejc (ref. 19) has shown reductions of approximately 
50 percent in clearance loss on axial turbine stages through the 
use of circumferential grooves in the rotor shroud. Although 
this concept does not appear to have been applied to a radial 
turbine configuration, the clearance phenomena are basically the 
same. The shroud and backface seal surface treatment is consid- 
ered viable for this application. Based on a 35-percent shroud 
and backface clearance loss reductions, a 1.35-point increase in 
performance is considered feasible. The total improvement pre- 
dicted for both stator loss reduction and rotor clearance 
effects is 2.35 points for the 1983-1985 time frame. 

Performance levels of the turbine inlet scroll, exhaust 
diffuser and deswirl vanes will be maintained at current 1979 
predicted levels. This does not imply that the performance 
potential of these components cannot be improved, but except for 
the deswirl vane, less research effort has been devoted to these 
components. In addition, the configurations ultimately adapted 
for this size shaft engine may differ from the representative 
configurations used in the parametric study. 

3.6 Preliminary Stator and Rotor Geometry 

The radial turbine stator preliminary geometry was estab- 
lished with the one-dimensional performance model. For a given 
stator trailing edge thickness of 0.030 in. (from ceramic con- 
siderations) and a trailing edge blockage of 10.0 percent (which 
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is consistent with the geometries used in the specific speed 
correlation) , the stator vane number and radial chord were 
established from a radial nozzle loading coefficient. The load- 
ing coefficient was based on the Zweiffel loading concept 
applied specifically to the radial nozzle geometry. 

The rotor inlet and exit radii and blade height were estab- 
lished with the one-dimensional performance model. The rotor 
flow path and axial length were evaluated from a meanline load- 
ing distribution through the rotor over a range of meridional 
flow shapes and axial lengths. The entire blade height through 
the rotor was treated as a stream tube and the continuity equa- 
tions, with an estimated rotor blade angle distribution, deter- 
mine the rotor mean-line relative velocity distribution. The 
condition of zero absolute vorticity was used to arrive at a 
blade-to-blade loading. 

3 . 7 Parametric Study Results 

A one-dimensional performance evaluation model for the tur- 
bine scroll, stage, and exhaust diffuser system was first estab- 
lished by modifying an existing radial turbine design program to 
account for the scroll and duct losses and projected 1985 per- 
formance improvements. The computer program calculates geo- 
metric and vector diagram quantities at the stator inlet, stator 
exit, rotor inlet and rotor exit mean-line stations. Due to the 
number of variables required to define the turbine flow path, 
the input to the program consisted of both free and fixed param- 
eters. The free parameters were: 

o Rotational speed (N) = 80,000, 95,000 and 110,000 rpm 

o Rotor inducer tip speed ^^t 3 ) = 549, 610, 701, and 782 

m/sec (1800, 2000, 2300, and 2565 fps) 

o Rotor inlet blade angle (/3„) = 0, 10, and 20 degrees 

D 

o Inducer tip-to-exducer tip radius ratio, 0.80 maximum 

2 

o Rotor exit annular area = 45.2 to 58.1 cm (7.0 to 9.0 
in. 2) 

o Rotor inducer work coefficient (which establishes an 
inducer-to-exducer work split, rotor exit swirl, and 
Mach number) 

o With and without deswirl vane 

o With and without predicted reaction effects 
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Fixed parameters used in the basic parametric study were: 

o Rotor blade number, set at 12 (The combined effects 
of rotor blade number and rotor exit blockage were 
examined for specific cases but are not included in 
the basic parametric curves) 

o Stator exit vane angle (defined at the rotor inlet 
radii) was set at 72 degrees 

o A stator loading parameter was fixed 

o The vane trailing edge thickness was fixed at 0.076 cm 
(0.030 in.) 

o Rotor clearances were set at 0.020, 0.020, and 0.010 
for rotor backface, axial, and radial, respectively 

o The vane trailing edge blockage was set at 10 percent, 
which was considered the maximum allowable for the 
specific speed correlation to be valid 

o The conditions of stator blockage, trailing edge 
thickness and loading criteria then set the stator 
vane number 

o The cycle conditions, except for efficiency and duct 
loss, were fixed 

Results of the basic parametric study are presented in Fig- 
ures 14 through 22. The basic characteristics of the aero- 
dynamic study are summarized below: 

o For a specified tip speed, the parametric study shows 
that an optimum system efficiency is obtained for all 
rotational speeds, both types of blading, with and 
without deswirl vanes 

o Rotor reaction significantly affects both the effi- 
ciency level and optimum inlet work coefficient at 
lower tip speed levels 

o Peak efficiency increases uniformly with increasing 
rotor inlet blade angle except at high tip speeds when 
rotor inlet incidence is eliminated 

o The level of tip speed required to achieve the system 
efficiency goal is relatively high [579 to 701 m/sec 
(1900 to 2300 fps)] 

o The turbine efficiency level increases with rota- 
tional speed; however, when the effect of rotor exit 
hub blockage is considered, this effect is reduced 
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SYSTEM TOTAL-TO-TOTAL EFFICIENCY, rjj_j 


• N = 80,000 RPM 

. % = 10 ° 

• NO DESWIRL VANE 

. Nb= 12 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 

ROTOR INLET WORK COEFFICIENT - X3 


Figure 15. Turbine System Parametric Study, N = 80,000 rpm, 

= 100. 
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SYSTEM TOTAL-TO-TOTAL EFFICIENCY, t? 


• N = 80,000 RPM 

• % = 20° 

• NO DESWIRL VANE 

• Nb = 12 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 


ROTOR INLET WORK COEFFICIENT, X3 


Figure 16. Turbine System Parametric Study, N = 80,000 rpm, 

= 20 °. 
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SYSTEM TOTAL-TO-TOTAL EFFICIENCY, j } j_j 


• N = 110,000 

• = 0 ° 

• Nb=12 



Figure 20. Turbine System Parametric Study, N = 110,000 rpm, jSg = 0°. 



SYSTEM TOTAL-TO-TOTAL EFFICIENCY, t]j_j 


• N = 110,000 RPM 

• = 10 ° 

• NO DESWIRL VANE 

• Nb = 12 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 


ROTOR INLET WORK COEFFICIENT, ^3 


Figure 21. Turbine System Parametric Study, N = 110,000 rpm, 

/3b = 10’. 
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SYSTEM TOTAL-TO-TOTAL EFFICIENCY, j 7 t_t 


• N = 110,000 RPM 

• = 20 ° 

• NO DESWIRL VANES 

• Nb=12 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 


ROTOR INLET WORK COEFFICIENT - X 3 


Figure 22. 


Turbine System Parametric Study, N = 110,000 rpm, 

/3g = 20". 
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o At lower rotor tip speeds [549 to 610 m/sec (1800 to 
2000 ft/sec) ] , system efficiency is increased from 2 
to 6 points with the addition of deswirl vanes 

In Figures 14 through 22, the peak efficiency for each 
value of inlet work coefficient analyzed is based on the maximum 
attainable efficiency for the range of rotor exit areas and 
radius ratios examined. Rotor exit areas of 45.2, 51.6 and 58.1 
cm^ (7, 8, and 9 in.^) were specified in the parametric study 
and, in general, peak efficiency occurred at or near the maximum 
value specified. The values of rotor exit area selected were 
established from preliminary rotor exit (zero swirl) continuity 
calculations. However, as shown in Figures 20 through 22, cer- 
tain combinations of rotation and tip speed will result in 
either negative or unacceptably low exit hub radii and no solu- 
tions will exist. Although high rotor exit area is desirable 
from a reduced exducer radial clearance loss standpoint, the 
higher exit area also requires increased hub blade thickness for 
a given blade number and blade stress level. The resultant 
increase in hub blockage could therefore negate the gains in 
clearance loss. 

For the radius ratio, peak system efficiency occurred at 
0.80 (the maximum allowable) for tip speeds up to 610 m/sec 
(2000 fps) . However, as optimum tip speed is approached, peak 
system efficiency occurs at lower radius ratios (0.65). These 
trends are associated with the rotor exit mean work coefficient. 
As tip speed is reduced from the optimum value, the level of 
rotor exit work coefficient increases rapidly. Therefore, the 
rotor exit conditions and duct loss are more sensitive to 
changes in the rotor exit mean radius. Under these conditions, 
the reduction in duct loss at higher radius ratios is greater 
than the increase in loss due to increased exducer clearance 
effects. It should be noted that the optimum radius ratio of 
0.80 for the lower tip speed cases also results in the highest 
rotor shroud curvature, and the performance effects for this 
parameter are currently not known. Therefore, if detailed mech- 
anical analysis indicates this level of tip speed is required to 
satisfy the rotor life goals, then a rotor internal flow solu- 
tion at both 0.80 and at lower values must be compared and the 
impact of the shroud loading evaluated. Figure 23 illustrates 
the effects of rotor exit area and radius ratio on peak system 
efficiency for 95,000 rpm and tip speeds of 610 m/sec (2000 fps) 
and 701 m/sec (2300 fps) . 

Specific vector diagrams and system parameters along the 
system efficiency envelope are presented in Figures 24 and 25 
for 80,000 rpm and the range of tip speeds investigated. Of 
particular note is the magnitude of the rotor exit relative 
critical velocity ratio associated with the deswirl vane cases 
at peak system efficiency. These higher values result in 
increased rotor exit blockage losses. In addition, since the 
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Figure 44. 


Campbell Diagram for 95,000-rpm Case with 45.2 cm2 
(7 in. 2) Exit Area. 
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Figure 42. Campbell Diagram for 95,000-rpm Case with 
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Figure 38. Campbell Diagram for 80,000-rpm Design 
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Figure 37. Tangential Thickness and Maximum Principal Stresses. 
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Figure 36. Two-Dimensional Finite Element Model for ISOPDQ. 
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The basic flow paths depicted in Figure 31 were examined 
for mechanical considerations using the 2-D finite element pro- 
gram, ISOPDQ (ref. 20), illustrated in Figure 36. Each flow 
path configuration was analyzed by an iteration on blade thick- 
ness and wheel/disk backface contour to achieve 2-D maximum 
principal blade stresses of 13.8 kN/cm^ (20 ksi) . 

Figure 37 illustrates the tangential blade thicknesses and 
maximum principal stresses for 80,000 rpm with 12 blades. This 
configuration exceeds the 24.1 kN/cm2 (35 ksi) stress limit in 
the disk. However, the blade vibration characteristics (Fig- 
ure 38) are acceptable. The most critical blade vibration con- 
dition expected for this engine application is the stator vane 
passing frequency. This size engine will be designed with a 
range of stator counts from 17 to 25. With a stator count of 23, 
this design has only one blade-natural frequency in the engine 
operating range. This possible frequency interference is an 
exducer torsional third mode (Figures 39 and 40) , which is not a 
problem with any existing radial turbines. 

Figure 41 illustrates the blade tangential thicknesses and 
the maximum principal stresses for 95,000 rpm, with 58.1 cm^ 
(9 in. 2) of rotor exit area and 12 blades. Both blade and disk 
stresses are within established limits. Based on the Campbell 
diagram of Figure 42, a stator count selection between 17 and 25 
would not reduce the number of blade natural frequencies below 
two in the operating range. To alter the blade for this speed to 
a shape similar to the 80,000-rpm case, which exhibits better 
vibration characteristics, an exit area of 45.2 cm2 (7 in. 2) was 
examined. Figure 43 illustrates the blade tangential thick- 
nesses and maximum principal stresses for this exit area. Fig- 
ure 44 shows that 23 stator vanes will reduce the number of 
blade-natural frequencies in the operating range to one. Fig- 
ures 45 and 46 illustrate the mode shapes for this blade config- 
uration. 

Figure 47 illustrates the blade tangential stresses, and 
Figure 48 depicts the blade-natural frequencies for 110,000 rpm. 
The maximum principal stresses in the disk are somewhat lower 
than for 95,000 rpm, but the blade vibration characteristics 
appear to add considerable risk for this wheel configuration. 
Although maximum principal disk stresses continue to decrease as 
physical speed increases (Figure 49) , blade vibration charac- 
teristics set a practical upward limit on physical speed to the 
95,000- to 100,000-rpm range. 

After the physical speed range was established, other mech- 
anical parameters were examined. The first of these variables 
was the number of blades. When the blade count was varied from 
10 to 16, very little change was observed in maximum stress at 
the center of the wheel (Figure 50) . These slight changes in 
bore stress with changes in blade count indicate that bore 
stress is more dependent on disk size than blade configuration. 
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with an axial blade length of 4.1 cm (1.6 in.). However, the 
final axial length and blade number was established during the 
detail rotor design. 

3 . 9 Parametric Study for Two-Dimensional Mechanical Analysis 

The mechanical analysis was initiated with the character- 
istic flow path configurations that resulted from the parametric 
study for 80,000, 95,000, and 110,000 rpm. 

The life objective of the program was to achieve a fast 
fracture modulus of rupture (MOR) probability of success of 
0.9999, using a selected ceramic material with a characteristic 
strength of 62.1 kN/cm^ (90 ksi) and a Weibull modulus of 15. 
These material property values are projections of selected 
ceramic material capability for the 1983 time period and were 
derived from a review of current ceramic matrials and associated 
processing capabilities. The following projected 1983 proper- 
ties for near-theoretical density Si 3 N 4 and SiC were used as a 
basis for this design study: 



Projected Strength (KSI) 

Weibull m 

Material 

RT 

1100“C 

RT 1100°C 

(2012°P) 

SI 3 N 4 

100 

90 

15 15 

Sic 

90 

90 

15 15 


To provide meaningful probability of success calculations, 
stresses must be integrated over both the surface area and vol- 
ume. Three-dimensional stress analysis is required to make the 
integrations. Initially, however, to examine a large number of 
configurations in the parametric study, a 2-D stress analysis at 
room temperature was used. 

Based on previous experience with ceramic radial turbines 
(utilizing similar material properties) , stress limits of 13.8 
kN/cm^ (20 ksi) in the blades and 24.1 kN/cm^ (35 ksi) in the 
disk were established. The blade limit for the 2-D analysis 
[13.8 kN/cm^ (20 ksi)] was fixed at a level providing adequate 
stress margin to accommodate stress concentrations in the com- 
plex geometry of the blade/disk interface region. The 2-D fin- 
ite element program does not compute concentrations in this 
transition region. 
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Figure 33. Rotor Solidity Study, Axial Length = 4.8 cm (1.9 in.). 
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3 . 8 Stator and Rotor Solidity Studies 


Since rotor inlet incidence loss is a function of the inlet 
stage work coefficient, the effect of stator exit flow angle can 
be examined independently at a constant work coefficient. The 
parametric study was based on a constant exit angle of 
72 degrees, which has proven to be a reasonable compromise 
between the sensitivity of manufacturing tolerances and rotor 
shroud and backface clearance effects. The effect of varying 
the stator exit angle is illustrated in Figure 32. For the 
study, stator trailing edge thickness and blockage were held at 
0.030 in. and 10 percent, respectively. The following charac- 
teristics are observed: 

o From continuity, increasing the stator exit angle 

increases the stator exit "b" width (also the rotor 
inlet ”b" width) 

o For a constant stator loading, exit blockage, and 

trailing edge thickness, increasing stator exit angle 
increases the stator aspect ratio (b/AR) 

o Increasing the stator exit angle reduces the rotor 

inlet relative velocity and increases rotor reaction 

o For an increase in stator exit angle from 66 to 

78 degrees, system efficiency increases by 1.3 points 
due to reduced shroud and backface clearance effects 
and increased reaction effects 

The rotor axial chord was evaluated on the basis of a mean 
line loading. The loading was established by treating the 
entire flow path as a stream tube and by solving the equations 
of continuity and energy from inlet to exit with the estimated 
flow path and rotor blade angle distribution. The condition of 
zero absolute vorticity was used to establish the blade suction 
and pressure surface velocities. 

The results of this analysis for the 701-m/sec (2300-fps) , 
95,000-rpm configuration are presented in Figures 33 to 35 for 
axial lengths of 4.8, 4.3, and 3.8 cm (1.9, 1.7, and 1.5 in.), 
respectively. For each axial length examined, a range of rotor 
blade numbers from 12 to 16 was evaluated. The results show 
that inducer loading was not affected by axial length, since the 
blade angle distribution in this region was fixed with radial 
blades. However, in the exducer and trailing edge region, the 
blade loading was significantly affected by axial length due to 
the rate of change in rotor blade angle. Conversely, the rotor 
blade number resulted in a change in rotor blade loading 
throughout the entire flow path, as expected. On the basis of 
this study, an axial length between 4.1 and 4.3 cm (1.6 and 1.7 
in.) resulted in satisfactory rotor blade velocity distri- 
butions. The majority of the mechanical analysis was performed 
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Figure 31. Characteristic Rotor Flow Path Configurations for 
Selected Parametric Study Rotational Speeds. 
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Figure 27. Peak System Efficiency for N = 95,000 rpm. 
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Figure 26. Peak System Efficiency for N = 80,000 rpm. 
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characteristic shown is at the rotor exit mean-line, signif- 
icantly higher relative critical velocity ratios would occur in 
the rotor exit tip region and shock losses should be evaluated. 

To define a tip speed envelope that will achieve the cycle 
performance goals, the peak system performance from Figures 14 
through 22 was first plotted as a function of tip speed and 
range of inducer blade angles, as shown in Figures 26 through 29. 
The results showed the magnitude of tip speed required for each 
rotational speed and also the effect of rotor inlet blade angle 
on required tip speed. Figure 29 shows these characteristics 
with deswirl vanes for 95,000 rpm. These results were 
cross-plotted in Figure 30 in terms of tip speed versus shaft 
speed, and a tip speed envelope was obtained. The envelope 
shows that without rotor exit blockage effects, tip speed 
requirements range between 579 and 716 m/sec (1900 and 2350 fps) 
depending on inlet blade angle, shaft speed and deswirl vane 
specification. 

Characteristic rotor flow paths for the three shaft speeds 
are presented in Figure 31. At 80,000 rpm, the basic flow path 
resulted in a relatively small blade, large disk configuration. 
Conversely, at 110,000 rpm the flow path resulted in a rela- 
tively large blade, small disk configuration. Although the 
radius ratio at 80,000 rpm could be decreased to match the exit 
hub radius at 110,000 rpm, this was not shown for two reasons. 
First, the parametric study showed increasing performance as 
radius ratio increased. Second, the inducer blade height 
increased by a factor of approximately 2.0, requiring a signif- 
icant increase in blade thickness to maintain equivalent inducer 
blade stress levels. In addition, the increased blade height 
was expected to result in undesirable blade vibration charac- 
teristics . 

Additional observations from Figure 31 are: 

o In addition to the smaller blade height dimensions 
associated with 80,000 rpm, the large disk also would 
increase rotor weight and inertia, which penalizes 
vehicle response time; therefore, driveability may 
suffer. However, rotor exit hub blockage would be 
minimized, and the number of blades could be 
increased . 

o At 110,000 rpm, rotor weight and inertia would be min- 
imized, and rotor blade height would be maximized, but 
the reduced rotor exit hub radius would significantly 
increase exit blockage or severely limit the number of 
blades. Mechanically, the relatively tall exducer 
blades could result in unacceptable blade vibration 
characteristics . 


48 



ROTOR EXIT MEAN 

ROTOR EXIT MEAN SWIRL ROTOR EXIT MEAN ABSOLUTE RELATIVE CRITICAL 

ANGLE,a, DEGREES CRITICAL VELOCITY RATIO, V/A. VELOCITY RATIO, W/A, 




ROTOR INLET WORK COEFFICIENT, Ag 


Figure 25. 


Turbine Vector Diagram and System Parameters 
Along Individual Tip Speed Envelopes, 

N = 80,000 rpm, = 0“ (Cont) . 


47 










SYSTEM TOTAL-TO-TOTAL EFFICIENCY, SYSTEM TOTAL-TO-TOTAL EFFICIENCY, 

•^-T/o^ ’^-T/0.6 



0.65 0.70 0.75 0.80 


RX 4 /RT 3 - EXDUCER TIP-TO-IIMDUCER TIP RADIUS RATIO 
(A) U ^3 = 610 m/sec (2000 fps) 



0.65 0.70 0.75 0.80 


Rt-/Rt« - EXDUCER TIP-TO-INDUCER TIP RADIUS RATIO 

^ w 

(B) Utg = 701 m/sec (2300 fps) 


Figure 23. Effects of Rotor Exit Area and Radius Ratio on Turbine 
System Efficiency, N = 95,000 rpra, / 3 „ = 0°, No Deswirl 
Vanes. 
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Figure 47. Tangential Thicknesses and Maximum Principal Stresses. 








Another parameter to be varied was tip speed. This vari- 
ation was accomplished by keeping the flowpath constant for 
95,000 rpm, with an exit area of 45.2 cm^ (7 in^.) and varying 
the inducer tip radius. Figures 51 and 52 show the tangential 
thicknesses and maximum principal stresses for tip speeds of 634 
m/sec (2080 fps) , and 70 m/sec (2300 fps) , respectively. The 
maximum disk stresses for the 701 m/sec (2300 fps) are somewhat 
higher than for the 634 m/sec (2080 fps) but are within the 
24.132 kN/cm^ (35 ksi) stress limit. Figure 53 summarizes the 
maximum disk stresses as tip speed is varied. 

3 . 10 Parametric Study of Detail Aerodynamic Analysis 

Since mechanical analysis indicated that disk stress and 
blade thickness would increase with the higher rotor exit area, 
and the rotor blade vibration characteristics were more favor- 
able at the lower 45.2 cm^ (7 in.^) rotor exit area, typical 
vector diagrams with and without deswirl vanes were selected for 
45.2 cm2 (7 in. 2) cases. Two-dimensional vector diagrams were 
calculated as a function of rotor exit blockage. 

For the case without deswirl vanes, the rotor exit average 
relative critical Mach number (W/Ac^) varies between 0.70 and 
0.80, and with a deswirl vane, the critical Mach number varies 
between 1.0 and 1.1. The general results in terms of a rotor 
relative total pressure loss ratio and hub blockage are pre- 
sented in Figure 54. When the effects of rotor blockage and 
blade number are combined, a minimum loss with and without 
deswirl vanes is defined as in Figure 55. The high penalty 
associated with the deswirl vane case is attributed to the 
higher rotor exit Mach number and exit angle due to the 
increased rotor exit swirl. 

To illustrate the combined effects of stage reaction and 
rotor exit blockage on system efficiency, the 95,000 rpm base 
solutions from Figure 17 were utilized. The results are 
presented in Figure 56. Since the parametric study indicated 
system efficiency increased with increasing rotor exit area, the 
overall effect of exit area, with reaction and blockage also was 
included. The radius ratio was fixed at 0.80. The charac- 
teristics derived from the basic parametric study are shown at 
the bottom of Figure 56 as a function of inducer tip speed and 
exit area. The change in efficiency as a function of stage 
reaction alone is next shown (derived from Figure 13) . As indi- 
cated, the benefits of higher exit area are already partially 
offset due to lower reaction. In addition, the higher incidence 
loss associated with low tip speeds is further compounded by 
relatively high losses due to reduced reaction. 

The rotor exducer blade thickness and blockage correlation 
derived from the 2-D mechanical analyses, together with the cal- 
culated 2-D rotor exit vector diagrams then were used to obtain 
a efficiency penalty for rotor exit blockage. Since the rotor 
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Figure 52. Tangential Thicknesses and Maximum Principal 
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Figure 54. Effect of Rotor Exit Hub Blockage on Rotor 
Loss. 






N = 9! 


/^B ~ ® 

Nb = 1 
RT4/RT3 = 0 


^FINAL 


SYSTEM EFFICIENCY 


S! 0.84 


I I 

"A4= 45.2 cm2 (7 |n.2) 
A4= 51.6 cm2 (8 in.2) 
Aa = 58.1 cm2 (g 2) 


+0.76 


ju +0.005 r^^BLOCKAGE SYSTEM 

U I I I 


< 

^ 0.00 
u 

3 -0.005 

CQ 

< - 0.010 


p50% BLOCKAGE 


-0.020 I REACTION 


g- 0.01 5 

H 

^- 0.010 

LU 

^ -0.005 
<1 

0.000 



<0 O.QOpSYSTEM EFFICIENCY, tjj_j 


Ai] = 58.1 cm 


I 

2 p 0.86 
ui . 

w ^ 

“I 


A4= 45.2 cm2 (7 j^ 2) 
A4= 51.6 cm2 (3 jp_2j 
A4= 58.1 cm2 (g jp_2j 


1800 1900 2000 2100 


cm/sec 


INDUCER TIP SPEED 


Figure 56. Effects of Reaction, Rotor Exit Area, and 
Blockage on System Efficiency. 







radius ratio was set at 0.80, the higher tip speeds result in a 
blockage penalty that is less than the reference turbine (i.e., 
positive Ai 7 blockage) . The combined effects of reaction and 
blockage are indicated by tj final. The advantages of the higher 
rotor exit area designs are reduced significantly and are lim- 
ited to higher tip speeds due to the maximum allowable blockage 
of 50 percent. 

Up to this point, the aerodynamic and mechanical analyses 
indicate the following trends; 

o The effects of reaction, exit area and blockage 
increased the required tip speed of all solutions by 
37 to 76 m/sec (120 to 250 fps) 

o The effect of rotor exit hub blockage increased the 
tip speed level of the 95,000-rpm solutions to the 
level required at 80,000 rpm. Since the disk stress, 
weight, and inertia are higher at 80,000 rpm, these 
solutions were eliminated 

o Comparison between the solutions at 95,000 and 110,000 
rpm show that disk stress is slightly lower at 110,000 
rpm, but the rotor exit hub blockage is significantly 
higher. To reduce the thickness requirement would 
require a reduction in rotor exit annular area to 
approximately 36.1 cm2 (5.6 in. 2), which decreases 
stage efficiency and increases the required tip speed. 
For this reason, the optimum solutions would be 
between 90,000 and 100,000 rpm. 

The tip speed characteristics for 95,000 rpm are presented 
in Figure 57 as a function of imposed effects of reaction, 
blockage, and exit annular area. The relative merit between 
radial and nonradial blading was examined next with a three- 
dimensional stress analysis. The rotor blade angle distribu- 
tions for a 10- and 20-degree inlet angle were estimated from 
previous nonradial rotor optimization. 
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Figure 57. Effect of Reaction, Blockage, and Exit Annular Area on Tip Speed 
Characteristics for Design System Efficiency of 0.87. 




4.0 DETAIL ROTOR AERODYNAMIC AND MECHANICAL DESIGN 


The three rotor configurations (without deswirl vanes) cor- 
responding to inlet blade angles of 0, 10, and 20 degrees were 
selected for detail aerodynamic and mechanical analysis. The 
tip speeds for these designs were based on the data presented in 
Figure 57 and represent the predicted values required to achieve 
the cycle system efficiency goal of 0.87 with all known losses 
and mechanical considerations accounted for. The design param- 
eters for the detail design configurations are presented in 
Table IV. Rotor flow path and blade geometries are presented in 
Appendix C. 

4 . 1 Aerodynamic Detail Design 

The rotor inlet and exit flow path dimensions (derived from 
Figure 57) were established from the optimized one-dimensional 
vector diagrams. The rotor solidity study, described in Section 
3.8, indicated that an axial length of 4.06 cm (1.60 in.) should 
result in satisfactory rotor blade loadings for the radial blade 
case. For the nonradial rotors, previous designs have demon- 
strated that additional axial length is required, due to the 
increased rotor turning. 

The detail rotor design was initiated with preliminary 
estimates of hub and shroud contours, blade thicknesses, and 
angle distributions. These quantities were input to a radial 
turbine geometry program along specified station lines (quasi- 
orthogonals) for the internal flow solution. The resultant data 
matrix was curve-fit to define the total blade geometry. The 
computer program is a modification of the current centrifugal 
impeller geometry program that allowed an arbitrary blade to be 
defined. 

To perform the internal flow analysis, the level and dis- 
tribution of losses in the rotor must be specified in addition 
to the rotor blade geometry and thermodynamic conditions. The 
magnitude of the total stage loss is determined from the effi- 
ciency analysis described in paragraph 3.4. The stator loss is 
estimated from stator reaction rig tests performed at 
AiResearch; the remaining loss is assigned to the rotor. 

The computer program for the rotor internal flow analysis 
solves the radial-equilibrium equation by satisfying the conti- 
nuity, momentum, and energy equations in the meridional plane in 
a manner similar to references 21 and 22. Blade surface veloci- 
ties are computed from the local rate of change in moment of 
momentum, the condition of zero absolute vorticity, and linear 
velocities between suction and pressure surfaces. Stanitz 
(ref. 23) has shown that this method produces satisfactory 
results compared with relaxation solutions of the potential flow 
equation. Numerous iterations between the geometry program and 
internal flow analysis program are required to achieve satisfac- 
tory blade loading for each thickness distribution examined. 
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TABLE IV. SELECTED CONFIGURATIONS FOR DETAIL AERODYNAMIC 
AND MECHANICAL DESIGN 


Design Parameters 

Units 

Configurations 



^ = 20“ 

i 

Rotational Speed N, 

rpm 

95,000 

95,000 

95,000 

Inducer Tip Speed, U. , 

1 

cm/sec (fps) 

701 

(2300) 

686 

(2250) 

671 

(2200) 

Stage Work Coefficient, 
^STG 


1.033 

1.091 

1.142 

Rotor Inlet Work 
Coefficient, 


1.000 

1.030 

1.085 

Specific Speed, N 

S 

rpm 

(ft^/Vsec^/^) 

0.537 

(69.20) 

0.537 

(69.20) 

0.537 

(69.20) 

Stage Reaction, R 3 ,pQ 


0.541 

0.545 

0.511 

Reynold ' s Number , 
Re X 10"^ 


1.028 

1.050 

1.074 

Exducer Tip-to-Inducer 
Tip Radius Ratio, 

^t /'^t 


0.632 

0.646 

0.661 

Rotor Exit Annular 
Area, A^ 


45.2 

(7) 

45.2 

(7) 

45.2 

(7) 

Rotor Exit Absolute, 
Swirl, 

degrees 

-0.76 

-13.40 

-11.50 

Number of Rotor Blades, 


12.0 

12.0 

12.0 

Turbine Inlet Scroll 
Loss , AP/P 


0.0057 

0.0056 

0.0056 

Rotor Exit Diffuser 
Loss , AP/P 


0.0303 

0.0322 

0.0311 

Stage Total-to-Total 
Efficiency, ^ 

^ ^ 1-4 


0.889 

0.892 

0.893 

j 

System Total-to-Total 
Efficiency, ^T-T Q_g 


0.870 

0.872 

































































The criteria for satisfactory blade loadings are monotoni- 
cally increasing velocities and minimum blade surface diffu- 
sions. For highly loaded rotors (that normally imply increased 
loading in the rotor inducer region) , these criteria are not 
satisfied. In addition, inlet flow deviations based on 
Stanitz's criteria are not expected to apply; therefore, the 
loadings in the inducer region become ill-defined. Under these 
conditions, more rigorous stream-function methods for evalu- 
ating the blade surface velocities are not warranted. The rotor 
internal flow analysis is further compounded by rotor exit flow 
deviations. Deviations are expected even with isentropic flows 
due to the rotor relative vorticity. Approximate methods for 
predicting these deviations were proposed by Mizumachi (ref. 24) 
in 1971, then Khalil (ref. 25) derived a rigorous method in 
1977. However, these deviations have been modified signifi- 
cantly by the transport of secondary flows from the inducer to 
the exducer tip region. 

The final rotor flow path configurations are presented in 
Figure 58. As indicated, the hub and shroud contours were main- 
tained for all three rotor inlet blade angle solutions. A con- 
stant flow path configuration allowed a more meaningful mechani- 
cal comparison, and previous nonradial rotor designs have pre- 
sented no particular advantage from modifying the rotor meridi- 
onal flow path under these conditions. Conversely, rotor axial 
length was fairly sensitive (to achieve satisfactory loading) 
when a rotor inlet blade angle was introduced. 

The change shown in Figure 58 is based on iterating axial 
length until approximately equal rotor inducer loadings for all 
three blade angles were achieved. An alternate approach would 
be to evaluate a rotor loading parameter along the hub, mean on 
shroud lines and iterate axial length until equal loading param- 
eters were achieved. The final rotor blade angle distributions 
for inlet blade angles of 0, 10, and 20 degrees are presented in 
Figure 59 as a function of percent meridional distance. Also 
shown are the locations at inlet and exit where flow deviations 
from the blade were assumed. 

The rotor inlet deviations are apparently due to rotor 
relative vorticity effects, and the location where the flow 
starts to follow the blade was based on Stanitz (ref. 1) corre- 
lation and was applied to radial turbines by Katsanis (ref. 22). 
The deviations at the rotor exit are even more complex, since 
the rotor relative vorticity effects are compounded by large 
secondary flow migration to the rotor exit tip region. The 
approach currently utilized at AiResearch is to assume the rotor 
exit flow is turned to the rotor mid-throat angle, and the addi- 
tional uncovered downstream turning is the total deviation at 
each radial location. The net effect with this procedure is to 
assume deviation increases from approximately zero at the hub to 
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a relatively high value at the tip and accounts for the devia- 
tions expected in the tip region from the secondary flows. How- 
ever, the parametric study results are based on precise 
inducer-to-exducer work splits. Achieving the exact vector 
diagrams desired would require either a significant increase in 
rotor internal flow knowledge or experimental iterations of 
rotor exit blade angles. 

The rotor suction and pressure surface velocity distribu- 
tions for the three inlet blade angle cases are presented in 
Figure 60 as a function of meridional distance. With 12 rotor 
blades and the adjustments in axial length, the loading through- 
out the blades is similar and is expected to achieve equal per- 
formance levels. 

4 . 2 Detailed Mechanical Design 

Ceramic materials possess large variations in strength 
properties, due to brittle characteristics. Consequently, prob- 
abilistic techniques must be utilized to accurately evaluate 
ceramic components. 

The probabilistic component evaluation technique employed 
by AiResearch is based on the Weibull statistical strength model 
(ref. 26). The working relationships permit the calculation to 
be made from a combination of local stress, component size, and 
Weibull parameters for both area- and volume-controlled frac- 
ture . 


CPF = 1 - e 


-R 


where R is defined as: 


and 


■ ■ • 4 /(- 


^oa\ 


Ba 


'ua 


V 


da 


A^ = Test-bar surface area 

= Test-bar volume 

= Elemental volumetric stress 

= Elemental surface stress 

0 = Volumetric zero failure stress 

ov 

0 = Surface zero failure stress 

oa 

= Volumetric characteristic stress 
= Surface characteristic stress 


Ba = Surface Weibull slope 
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Bv = Volumetric Weibull slope 
da = Differential surface area 
dv = Differential volume 

This probabilistic computer code distinguishes between surface 
and volume strength relationships and permits assessment of 
overall component risk in a method compatible with finite 
element techniques. All three principal stresses are considered 
and the codes can assess effect of multiaxial stress states and 
change in reliability due to proof testing (ref, 27). 

As previously mentioned, cumulative probabilities of suc- 
cess computations require 3-D finite element stress results. 
The 3-D model for the optimized radial-blade configuration is 
shown in Figure 61. Detailed stress analyses for a similar 
radial inflow turbine rotor configuration in the AFTlOl program, 
under DOE/NASA contract DEN3-167, have shown that stresses and 
cumulative probabilities of success are quite similar for both 
room- temperature and maximum power steady-state operating con- 
ditions. This AGT study also showed that minimum cumulative 
probabilities of success are computed during transient opera- 
tion. Unfortunately, transient thermal studies could not be 
completed for the radial and non-radial blade configurations of 
this study because of financial restrictions. Consequently, all 
comparative stress analyses were carried out for room- 
temperature operation. However, probabilities of success were 
computed using anticipated 1983 material properties at elevated 
temperature (characteristic strength of 62.053 KN/cm2 (90.0 KSI) 
and a Weibull modulus of 15.0) . The parametric study indicated 
that the 95,000-rpm wheel, with an inlet tip speed of 701 m/sec 
(2300 fps) , was the preferred radial-blade configuration and is 
summarized on Table V. When this configuration was first exam- 
ined mechanically (Figure 52) , disk stresses were marginally 
high. To reduce disk stresses, the saddle height of the wheel 
was reduced by removing disk material above the self-sustaining 
radius of the disk. The resulting wheel stress field at room 
temperature is illustrated in Figure 62. The peak bore stress 
was reduced from 24.132 kN/cm2 (35.0 ksi) to 22.339 kN/cm^ 
(32.4 ksi). The resulting probability of survival for the rad- 
ial wheel using anticipated 1983 ceramic material properties 
[characteristic strength of 62.053 kN/cm^ (90.0 ksi) and 
Weibull modulus of 15.0] is 0.9999. This cumulative probability 
of survival (CPS) meets the objective of the program. 

Once the baseline radial blade design was analyzed, the 
nonradial inlet design feature was studied using the 3-D stress 
and probabilistic evaluation computer code. Initially, an inlet 
angle of 10 degrees was examined with the blade thickness dis- 
tribution for the radial blade configuration. To achieve the 
same efficiency as the radial blade design, the required tip 
speed for this nonradial configuration was reduced to 686 m/sec 
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Figure 61. Typical Three-Dimensional Finite Element 
Model. 
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TABLE V. OPTIMUM CONFIGURATION DERIVED FROM PARAMETRIC 
STUDY AND 3-D MECHANICAL ANALYSIS 


Parameter 

Optimum 

Shaft speed 

95,000 

Inducer tip speed 

701 cm/sec 
(2300 fps) 

Rotor inlet blade angle 

0 degrees 

Rotor blade number 

12 to 14 

Rotor exit swirl 

-10 degrees 

Rotor exit critical velocity ratio 

0.449 

Rotor inducer-to-exducer tip 
radius ratio 

0.632 
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(2250 fps) . Principal stress results for this configuration 
(Figure 63) indicate high bending stresses in the wheel inlet. 
These produce a low probability of survival of 0.9394. Experi- 
ence with nonradial blades from other research programs indi- 
cated that these bending stresses could be reduced by incorpo- 
rating rake to the blade design. Figures 64 and 65 illustrate 
configurations without and with rake (-30 degrees ) , respec- 
tively. The principal stresses for the rake case (Figure 66) 
improve the cumulative probability of success to 0.9986. 
Although this CPS value is significantly better than the no-rake 
design, the goal of 0.9999 was not reached. 

Blade thickness distribution was the next parameter to be 
studied. Increased blade thickness was incorporated in the 
inducer portion of the -30 degree rake case (Figure 67) . The 
resulting principal stresses, illustrated in Figure 68, produce 
a cumulative probability of survival of 0.9993. 

Based on the improvements in blade stresses and CPS with 
the modified blade thickness distribution, a new distribution 
was generated that increased the blade thickness in both the 
inducer and exducer. The resulting principal stresses (Fig- 
ure 69) created a cumulative survival probability of 0.9998. 
The disk stresses were increased slightly due to thicker blades. 

The thickness distribution was modified once again by 
thickening the blades in the high stress regions. The resulting 
blade stresses were significantly reduced (Figure 70) . Disk 
stresses were again increased slightly. The survival rate for 
this configuration, 0.9997, was lower than the previous blade 
thickness design because of the higher disk stresses. 

The 20-degree inlet configuration was analyzed with two 
different blade thicknesses: the radial blade and the thickest 
distribution for the 10-degree inlet. To achieve an efficiency 
equivalent to the radial design, tip speed was reduced to 671 
m/sec (2200 fps). As expected, the radial blade thickness dis- 
tribution produced extremely high blade stresses (Figure 71) and 
an unacceptable cumulative probability of survival of 0.0001. 
With the thick-blade distribution, the blade stresses were con- 
siderably reduced (Figure 72) , and resulting cumulative prob- 
ability of success was 0.4865. The blade stresses could be 
reduced by a further increase in blade thickness, but these 
increases would increase disk stresses. Increased thickness 
would prevent the 20-degree inlet case from ever achieving a 
survival rate as high as the radial blade design. 

The peak blade and disk stresses and cumulative probabili- 
ties of survival for the different wheel geometries are tabu- 
lated on Table VI. In summary, the radial blade is the only wheel 
design that meets the 0.9999 CPS program goal; nonradial config- 
urations do not possess survival rates as high. Blade thick- 
nesses had to be increased significantly to keep the blade 
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SCALE FACTOR = 1.50 RPM = 95,000 
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Figure 63. Maximum Principal Stresses, 
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Figure 64. Blade Shape, Original Thickness (Rake 
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NONRADIAL. BETA = 10 
RAKE = -30 



Figure 65. Blade Shape, Original Thickness (Rake = -30). 



NON RADIAL. BETA = 1C 
RAKE = -30 


SUCTION SIDE 


64.5 KSI 


^^31 0 KSI 

ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 

SCALE FACTOR = 1.50 
<D MAXIMUM 4^ MINIMUM 


31 0 KSI 

ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 
RPM = 95,000 

TIP SPEED = 686 m/sec (2250 fps) 


Figure 66. Maximum Principal Stresses, Original Thickness. 
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Figure 67- Blade Shape, Modified Thickness 
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ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 


ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 


SCALE FACTOR = 1.50 RPM = 95,000 

<!> MAXIMUM A MINIMUM TIP SPEED = 686 m/sec (2250 fps) 


Figure 68. Maximum Principal Stresses, Modified Thickness (1) . 
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ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 

SCALE FACTOR = 1.50 RPM = 95,000 

O MAXIMUM A MINIMUM TIP SPEED = 686 m/sec (2250 fps) 


Figure 69. Maximum Principal Stresses, Modified Thickness (2). 









ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 

SCALE FACTOR = 1.50 RPM = 95,000 

O MAXIMUM A MINIMUM TIP SPEED = 686 m/sec (2250 fps) 


Figure 70. Maximum Principal Stresses, Modified Thickness (3). 
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ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) ISOPLETH INTERVAL = 3.447 kN/cm^ (5 ksi) 

SCALE FACTOR = 1.50 RPM = 95,000 

<!> MAXIMUM A MINIMUM TIP SPEED = 671 m/sec (2200 fps) 


Figure 72. Maximum Principal Stresses, Modified Thickness (3). 




TABLE VI. STRESS AND CUMULATIVE PROBVBILITY OF SUCCESS SUMMARY AT 95,000 RPM 


Configuration 

Rake 

Blade Inlet 
Tip Speed, 
ra/sec 
(fps) 

Max Blade 
Stress , 
kN/cm2 
(ksi) 

Max Disk 
Stress , 
kN/cm^ 
(ksi) 

Cumulative Probability 
of Success 

Contribution 

Total 

from Blades 

from Disk 

Radial Blade 



701 

16.8 

22.3 

1.00000 

0.99989 

0.99989 




(2300) 

(24.4) 

(32.4) 




10-Degree Inlet 

0 


686 

59.4 

21.4 

0.93950 

0.99992 

0.93942 

Radial Blade Thicknesses 



(2250) 

(86.1) 

(31.0) 





-30 



44.5 

21.4 

0.99862 

0.99992 

0.99855 






(64.5) 

(31.0) 




Modified Thickness 





33.8 

22.0 

0.99940 

0.99987 

0.99928 

Distribution #1 





(49.0) 

(31.9) 




#2 





31.4 

22.5 

0.99996 

0.99981 

0.99976 






(45.6) 

(32.7) 




#3 





28.1 

23.0 

0.99996 

0.99972 

0.99970 






(40.7) 

(33.3) 




20-Degree Inlet 



671 

91.5 

21.4 

0.00006 

0.99992 

0.00006 

Radial Blade Thickness 



(2200) 

(132.7) 

(31.1) 




Modified Thickness 



671 

52.4 

23.0 

0.71518 

0.99972 

0.71498 

Distribution #3 



(2200) 

(76.0) 

(33.3) 







stresses at acceptable levels. Increased blade load increased 
disk stresses and caused disk survival rates to be lower than 
those for the radial blade design. With sufficient blade thick- 
ness and rake adjustments, moderate amounts of backsweep (10 
degrees or less) can be tolerated with nearly the same cps as a 
radial bladed wheel. 

The 10-degree inlet design that possessed the highest cumu- 
lative probability of survival (modified thickness No. 2) and 
the radial blade configuration were compared, using ceramic 
properties more representative of present day materials [a char- 
acteristic strength of 44.816 kN/cm2 (65 ksi) and a Weibull mod- 
ulus of 12]. With these lower values, the radial wheel again 
possessed a better survival rate — 0.9686 versus 0.9359 for the 
nonradial design. 
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5.0 OFF-DESIGN PERFORMANCE EVALUATION 


For a single-shaft configuration, turbine operating condi- 
tions will vary significantly from maximum power to idle. 
Therefore, the turbine design derived from the maximum power 
optimization study may not provide minimum SFC and maximum 
driveability over the entire duty cycle. Although the turbine 
design point could be selected at any power setting and then 
examined off-design, this approach is not necessary. Evaluating 
the off-design characteristics at different "simulated" power 
settings can be accomplished merely by selecting a range of vec- 
tor diagrams at maximum power setting. For example, since tur- 
bine flow, speed, and pressure ratio will be reduced at idle, 
high positive swirl would be expected and, with the duct loss 
characteristics used in the parametric study, would result in 
decreased performance. Of course, turbine performance at idle 
is also dependent on a number of other parameters, such as work 
coefficient and reaction. 

The point is that if the vector diagram at maximum power is 
changed (for example, from -10.0 to -40.0 degrees exit swirl), 
all off-design vector diagrams will shift, and new efficiency 
characteristics will be defined. In addition, the system effi- 
ciency envelope already provides alternate vector diagrams that 
would maintain maximum efficiency at the design-point. 

The effect of design-point vector diagram selection on 
overall duty cycle performance will be illustrated based on the 
results of the current NASA-DOE Advanced Gas Turbine Power Train 
System Development Program (Contract No. DEN3-167) . Since the 
cycle and selected design point parameters are very similar to 
the 701 m/sec2 (2300 fps) tip speed radial-bladed configura- 
tions, the results are directly applicable to the present study. 

Figure 73 shows three selected inlet work coefficients (and 
the corresponding rotor exit swirl levels) on the system effi- 
ciency envelope for 95,000 rpm and a tip speed of 701 m/sec 
(2300 fps). The three cases are designated low (-9.8 degrees) 
swirl, high (-27.0 degrees) swirl, and maximum (-37.0 degrees), 
and the low swirl case corresponds to peak system efficiency at 
maximum power. 

Figure 74 shows the results of each vector diagram over the 
entire duty cycle. At idle, the maximum swirl design increased 
efficiency by approximately 6.0 points compared to the low-swirl 
design. At maximum throttle and cruise, the performance of the 
high-swirl and maximum-swirl designs was almost equivalent. At 
maximum power, the maximum swirl efficiency was reduced by 
approximately 2.5 points. Overall, the high-swirl design 
appeared to offer the best compromise between reduced maximum 
power and increased part-power characteristics. However, cycle 
studies 'will be required for each case before a final selection 
can be made. Not only were the off-design characteristics 
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Figure 73. Maximum Power Design -Point Selection for Off-Design 
Performance Evaluation. 
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Figure 74. Off-Design Performance Characteristics as a Function 
of Design-Point Exit Swirl (from NASA-DOE Advanced 
Gas Turbine Power-Train System Development Program) . 
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strongly dependent on the exhaust duct swirl loss charac- 
teristics (and, therefore, had to be known before a vector dia- 
gram was selected) , but the turbine Reynolds number was reduced 
significantly at part-power and had to be accounted for in the 
off-design analysis. 

The off-design analysis also confirms that the use of 
deswirl vanes is not viable for the single-shaft engine applica- 
tion due to the large rotor exit swirl excursion during the nor- 
mal duty cycle operation. The parametric study has shown that 
deswirl vanes can significantly reduce downstream duct losses; 
therefore, deswirl vanes should not be ignored during the tur- 
bine optimization process. 



6.0 CONCLUSIONS 


The results of this program show that exceptionally high 
radial turbine tip speeds are feasible with the projected 1983 
ceramic material properties (62.1 KN/cm^ characteristic 
strength and a Weibull modulus of 15) . These results indicate a 
high-work ceramic radial turbine with 0.076 cm (0.030 inch) 
thick blade tips could achieve both the cycle performance goal 
of 0.87 ^ ip-T rotor mechanical integrity goal of 0.9999 

cumulative probability of success for a single-shaft automotive 
engine application. 

Three-dimensional rotor analyses showed that viable high 
tip speed designs were derived from a design technique that com- 
bined the characteristic of low disk volume and relatively deep 
scallops with a blending of the blades and disk in the rotor 
back-face region. With radial blade element rotors, three- 
dimensional mechanical analysis showed that the cumulative prob- 
ability of success (CPS) goal of 0.9999 was achieved at a tip 
speed of 701 m/sec (2300 fps) . 

With a rotor inlet blade angle of 10 degrees, the rotor tip 
speed requirement was reduced to 686 m/sec (2250 fps) , and 
detailed mechanical optimization of this design showed that the 
CPS was reduced to 0.9998. At a tip speed of 671 m/sec 
(2200 fps) and a rotor inlet blade angle of 20 degrees, addi- 
tional blade thickness was required to achieve tolerable blade 
stress levels. The increased loading from this additional blade 
material caused higher disk stresses that adversely effected the 
CPS value. Therefore, the conclusion is that although nonradial 
rotor blading allowed a reduction in tip speed, this advantage 
was offset by an increase in rotor blade stress levels. This 
conclusion also would apply with reduced ceramic material prop- 
erties; that is, the unique ceramic characteristic of low duc- 
tility and density favors the lower global stress field asso- 
ciated with the radial bladed designs. Therefore, if the pro- 
jected ceramic properties are not achieved, turbine performance 
potential would be reduced with both radial and nonradial blad- 
ing, and fuel economy would suffer. 

Additional observations, resulting from the advanced 
radial turbine program, are: 

o The turbine system optimization technique more accu- 
rately identified peak system efficiency configura- 
tions compared to an analysis approach involving indi- 
vidual stage and exhaust duct optimization. 
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o 


The estimated effects of reaction and rotor exit bloc- 
kage were relatively important for a large number of 
solutions examined. Additional work is desirable to 
more accurately refine/define performance correla- 
tions for these effects. 

o For engine operation requiring only relatively small 
variations from turbine design-point operating condi- 
tions, the analysis showed that significant improve- 
ments in performance were provided with the applica- 
tion of deswirl vanes. However, for a single-shaft 
automotive engine, the large excursion in rotor exit 
swirl eliminated consideration of deswirl vanes for 
this application. 

o The off-design analysis showed that the integrated 
duty cycle performance was highly dependent on the 
selected design-point vector diagram. This implies 
that either the effect of rotor exit swirl on duct 
performance must be accurately known in advance of the 
study analysis, or that several design-point vector 
diagram rotor configurations must be evaluated exper- 
imentally to arrive at the optimum duty cycle config- 
uration. 
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APPENDIX A 


ABBREVIATIONS AND SYMBOLS 
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APPENDIX A 


a 

b 

b/AR 

c 

g 

q 

w 

Wv /^/6 

A 

J 

L 

M 

N 

R 

U 

V 

W 

Z 


ABBREVIATIONS AND SYMBOLS 
Scroll inlet diameter, cm (in.) 

Stator exit and rotor inlet passage width, cm (in.) 
Stator aspect ratio 
ln(r^ + a)/r^ 

Force-mass conversion constant, 32.174 Ibm-ft/ 
lbf-sec2 

2 2 

Dynamic pressure, kN/m (Ib/in. ) 

Scroll inside radius, cm (in.) 

Stage mass flow, kg/sec (Ib/sec) 

Stage inlet corrected flow, kg/sec (Ib/sec) 

2 2 

Annular area, cm (in. ) 

Mechanical equivalent of heat, 778.0 ft-lbf/Btu 

Rotor axial length, diffuser length, cm (in.) 

Meridional length, cm (in.) /Rotor exit mean, cm 
(in. ) 

Rotational or shaft speed, rpm 
Radius, cm (in.) 

Blade speed, m/sec (ft/sec) 

Absolute velocity, m/sec (ft/sec) 

Relative velocity, m/sec (ft/sec) 

Axial dimension, cm (in.) 
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ABBREVIATIONS AND SYMBOLS (Contd) 

1 2 

A^j. - Critical velocity gRgT,pl , m/sec (ft/sec) 
AR - Deswirl vane aspect ratio 


Ca - Axial rotor shroud clearance, cm (in.) 

Cb - Rotor backface clearance, cm (in.) 

Cr - Radial rotor shroud clearance, cm (in.) 

Djj - Hydraulic diameter, m (ft) 

Max - Maximum 

MC - Midchannel 

Min - Minimum loss coefficient 

M„ - Absolute flow Mach number 
N 

M„ - Relative flow Mach number 


N/# 


N 


B 


N, 


N. 


R, 


T-DE 

^R 
PS 


Corrected speed, rpm 
Rotor blade number 


Specific speed 






\l/2 




\3/4 


I SEN/ 

Number of stator vanes 


rpm (ft^/^/sec^/^) 


Total-to-dif fuser exit static pressure ratio 
Pressure ratio 


Pressure surface 

2 2 

Static pressure, kN/m (Ib/in. ) 


p 

STD 


P 


T 


2 

Standard atmospheric pressure, 1.01325 kN/m 
(14.696 lb/in.2) 

2 2 

Total pressure, kN/m (Ib/in. ) 
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ABBREVIATIONS AND SYMBOLS (Contd) 

Re - Reynolds number for turbine stage (w/R^ ) 

Re - Reynolds number for turbine scroll 

R - Gas constant, J/kg K (f t-lb/lb-®R) 
y 

^STG ” Stage reaction 

S.F. - Slip factor for radial blades 

SS - Suction surface 

*^STD ” Standard atmospheric temperature, 288.2 K (518. 7"R) 
TPLP - Total pressure loss parameter 


Rp - Diffuser recovery 


- P, 


T^ - Absolute total temperature, K (°R) 

- Blade tip speed, m/sec (ft/sec) 

- Absolute radial velocity, m/sec (ft/sec) 

- Absolute tangential velocity, m/sec (ft/sec) 

V - Absolute axial velocity, m/sec (ft/sec) 

W - Relative velocity, m/sec (ft/sec) 

W/A^j. - Relative critical velocity ratio 

OL - Rotor exit swirl angle, absolute flow angle, 
degrees (measured from axial direction) 

/3 - Rotor blade angle, degrees 

j3g - Rotor inlet blade angle, degrees 

y - Ratio of specific heats 
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ABBREVIATIONS AND SYMBOLS (Contd) 

6 - Ratio of turbine inlet total pressure to standard 
atmospheric pressure 

Aa - Total stator turning, degrees 

Ah - Static enthalpy change across rotor, kJ/kg (Btu/lb) 

AH^ - Specific stage work, kJ/kg (Btu/lb) (based on total 
conditions) 

AH^/0 - Specific stage corrected work, kJ/kg (Btu/lb) 
(based on total conditions) 

AP/P - Total pressure loss 


AR - Stator radial chord, rotor exit blade height, 
cm (in.) 

- Change in efficiency 
^ - Efficiency 

9- Ratio of turbine inlet total temperature to 
standard atmospheric temperature 

- Equivalent cone angle defined in Appendix B 
eq 

X - Work coefficient (V^/U) 

X - Friction resistance coefficient for scroll loss 

gJAH^/U^ 

Xo j - Ideal rotor inlet work coefficient 


- Stage work coefficient 



^3 , act 


- Actual rotor inlet work coefficient from vector 
diagram 


X 4 M " Rotor exit mean work coefficient 

2 2 

u - Kinematic viscosity, m /sec (ft /sec) 


2 2 

M - Absolute viscosity, N sec/m (Ib-sec/ft ) 


O) 


Loss coefficient 
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ABBREVIATIONS AND SYMBOLS (Contd) 

SUBSCRIPTS 

base - Base turbine total efficiency from specific speed 
curve 

blockage - Rotor exit blockage 
D - Diffusion efficiency 

i - Turbine stage total efficiency with incidence loss 

final - Total system efficiency with blockage and reaction 
effects 

ISEN - Isentropic 

loss - Rotor relative total pressure loss 
max - Maximum 
min - Minimum 
reaction - Stage reaction 

T-DE - Total to diffuser exit static 
T - Tip 

T-sys - Inlet total to exhaust diffuser exit total 
T-T - stage or system total-to-total efficiency 

Turbine System Station Nomenclature 

0 - Scroll inlet (combustor discharge) 

1 - Stator inlet 

2 - Stator exit 

3 - Rotor inlet 

4 - Rotor exit (deswirl vane inlet) 

5 - Deswirl vane exit (exhaust diffuser inlet) 

6 - Exhaust diffuser exit 
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APPENDIX B 


DE SWIRL VANE LOSS MODEL 
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APPENDIX B 


DESWIRL VANE LOSS MODEL 


For highly loaded radial turbine designs. Section 3.7 
showed that system performance is significantly increased with 
the application of rotor exit deswirl vanes. The concept is to 
relieve the rotor inducer loading by increasing rotor exit work 
coefficient, which in turn increases rotor exit velocity and 
swirl levels. The deswirl vane is then utilized to reestablish 
axial flow and reduce velocity by the difference between the 
absolute and axial components. 

A comprehensive evaluation of deswirl vane losses was con- 
ducted by Mitchell and Soileau (ref. 3) between 1974 and 1977 
under Air Force sponsorship. From a preliminary design stand- 
point, the major results of this program are; 

o Non-series - Airfoils demonstrated a reduction in 
total pressure, higher diffusion efficiencies, and 
reduced tendencies for flow separation relative to 
series airfoils. 

o Conical angle - An equivalent conical angle correla- 
tion (^eq) was extended for diffusing cascades into a 
range of cone angles applicable to advanced engine 
turbine exit guide vanes. The correlation is used to 
define optimum diffusion efficiencies and was suc- 
cessfully demonstrated in 3-D annular cascade test- 
ing. 

The correlation derived by Mitchell and Soileau (from dif- 
fusing cascade data) is presented in Figure 75. The correlation 
relates equivalent cone angle (the rate of diffusion) with dif- 
fusion efficiency {Vj)) as a function of equivalent area ratio 
(amount of diffusion) . 

The equivalent cone angle is defined with the nomenclature 
of this report as; 


S = tan 
^eq 


-1 


1/2 \ l/2 

0. 564^^1 ^cos - cos a^j 


With parallel end-walls, the equivalent area ratio is: 

^5^'^4 ~ cos a. 
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2 4 6 8 10 12 14 16 

EQUIVALENT CONE ANGLE - DEGREES 


Figure 75. Equivalent Conical Angle Correlation of 
Diffusing Cascade Data (ref. 3) . 


The diffusion efficiency relates the actual- to-ideal 
recovery in the following manner: 


{^S5 “ ^ 34 ) 

'/ 

C^T4 ^34 (measured). 

1 

1^35 ^34 ) 


,^T4 “ ^34 (ideal)] 


For the advanced radial turbine parametric study, the line 
of optimum diffuser efficiency shown in Figure 75 was used to 
define a relationship between equivalent area ratio and effi- 
ciency (Figure 76) . Since parallel end-walls Were assumed for 
the study (consistent with the selected diffuser geometry) , the 
rotor exit swirl angle is shown as a function of area ratio. 
From the data of Figure 76, rotor exit swirl levels from 37.47 
to 50 degrees can be evaluated. Fifty degrees appeared reason- 
able for the upper limit and, therefore ,• the maximum rotor exit 
angle considered with deswirl vanes was set at this value. How- 
ever, at the lower swirl levels, the trade-off occurrence between 
deswirl and no deswirl vanes is not clear. Therefore, the data 
from Figure 76 was extrapolated to 30 degrees, and then conven- 
tional axial compressor cascade correlations were used to define 
the deswirl vane loss to 0-degree swirl. 

The loss model for 0- to 30-degree exit swirls was based on 
an in-house correlation between a total loss parameter (TPLP) 
and diffusion factor (Df) with corrections for end-wall effects. 
With representative values of solidity (a) and diffusion factor 
of 1.4 and 0.4, respectively, a loss coefficient (co) was calcu- 
lated with the following relationship: 

- 2 (TPLP) 

(j) = - 

cos 


For exit swirl levels from 30 to 50 degrees, nonseries air- 
foils were^ assumed to apply, and the total pressure loss for the 
deswirl vanes were calculated. 


From the one-dimensional rotor vector exit mean-line vector 
diagram, the following quantities were known: 



“4 


^T4 ' ^S4 


and 


w/t 
A P,p 


4 


The deswirl vane ideal exit conditions can be obtained from 
the area ratio and continuity relationships. The equation is 
rearranged to solve for the deswirl vane exit static pressure 
(since diffusion efficiency is known) : 


P = , S5 ~ ^S 4 

^S5 (measured) = ^d\p ^4 - P 54 


ideal 


[^T4 ^S4 (measured)] 
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DIFFUSION EFFICIENCY, ’Jp ROTOR EXIT ABSOLUTE SWIRL ANGLE, a^. DEGREES 



EQUIVALENT AREA RATIO, A5/A41 

'eq 


Figure 76. Deswirl Equivalent Area Ratio for Maximxam 
Diffusion Efficiency (ref. 3). 
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The deswirl vane exit measured static pressure is used to 
solve for continuity, and an exit measured total pressure is 
defined . 

The downstream diffuser loss is based on the axial velocity 
(zero swirl), total pressure, and loss coefficient established 
at the deswirl vane exit. The summation of the deswirl vane and 
diffuser losses are used with the stage efficiency to define the 
system efficiency for the particular configuration being exam- 
ined in the parametric study. 
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APPENDIX C 


ROTOR FLOW PATH AND BLADE 
GEOMETRY FOR /3g = 0 , 10 AND 20“ 
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APPENDIX C 


ROTOR FLOW PATH AND BLADE 
GEOMETRY FOR = 0, 10 AND 20“ 


The rotor flow path and blade geometry for the three rotor 
inlet blade angle configurations are presented in Tables I, II 
and III for (3q - 0, 10 and 20“ , respectively. The data array 
included is for 7 equally-spaced modal points along a quasi- 
orthogonal station line. The first station point is on the 
shroud streamline and the last (seventh) point defines the hub 
streamline. The blade normal thickness is defined as the base- 
line distribution for = 0“. The variables for each station 
line are defined below: 

R - Local radius on station line (inches) 

Z - Local axial dimension on station line (inches) 

(Z = 0 is defined at rotor trailing edge) 

BETA - True blade angle (tan /3 = ^^) 
where; 

(f) = Polar angle about the axis of rotation 
m = Meridional distance along a streamline 

DEV - deviation; the difference between the rotor flow 
and blade angle in degrees. Rotor inlet and exit 
deviations were applied in the internal flow solu- 
tion 

THICK - Blade normal thickness (inches) 

The final thickness array used to minimize the blade bend- 
ing stresses for fio = 10 and 20“ is presented in Tables IV and V. 
Station data in tnese two tables starts at the blade inlet, 
instead of the exit, and contains tangential thickness in place 
of normal thickness. Beta distribution for these thickness 
distributions is similar to the distributions in Tables II and 
III. Dimensions are in inches. 
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Rot. Number 4 = Blade Trailing Edge 
Number 19 = Blade Leading Edge 


ST 

R 

I , 
BETA 1 
DEV .A 
.JHICK Ig 

! ST 
.R 

BETA 1 
DEV ,, 
jTHICK 

ST 


. 5<)300 

.70900 

.67000 

.00000 

.03803 


.65300 

.83600 

.18000 

.00000 

.03816 


'IR 

z 

BETA 
DEV . 
'.THICK g 

i ST 

R — - 

Z , 
BETA 1 
DEV 1 
THICK D 

! ST 

R 

BETA 2 
DEV \ 
THICK 4 

ST 

R 

Z , 

BETA-2 

DEV } 

■'jTHICK 1 

S“ 

R 

..beta 2 

DEV 1 
■THICK - 

1 D 

i ! S 

R 

Z , 
- BETA 2 
' DEV 1 
'^thicko 

S 

R 

2 T 

BETA-2 
DEV 1 
THICK® 

S 

R 

^‘l|TA-2 

DEV 1 
2.THICK0 
■ 0 


.73100 

.97800 

.43000 

.00000 

.04046 


.85000 

.11500 

.49000 

.00000 

.04758 


.01400 

.24300 

.20000 

00000 

.05690 


.19000 

.33300 

.12000 

.00000 

.05579 


.39300 

.39300 

.16300 

.00000 

.04252 


.59300 

.42300 

.00000 

.00000 

.03498 


.77430 

.43850 

.00000 

.00000 

.03013 


,92390 

.43850 

.00000 

.00000 

.03596 


1.50300 

1.41200 

1.32100 

1.23000 

.79100 

.87300 

.95500 

1.04000 

24.07000 

17.47000 

10.87000 

4.27000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.05136 

.07159 

.09816 

.15148 

1.57900 

1.50500 

1.43200 

1.35800 

.93000 

1 .02200 

1.11300 

1.20500 

18.38000 

12.58000 

6.78000 

.98500 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.05302 

.07266 

.10703 

.15625 

1.67700 

1.62200 

1.56800 

1.51300 

1.07400 

1.16900 

1.26500 

1.35780 

12.36000 

8.29000 

4.21000 

.14000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0. 00000 

.05711 

.08627 

.12351 

.17630 

1.81250 

1 .77500 

1.73750 

1.70000 

1.20300 

1.29000 

1.37800 

1.46500 

- 6 . 38000 - 

- 4.26000 

2.14000 

.02000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.06803 

.10223 

.13418 

.18345 

1.99300 

1.97200 

1 .95100 

1.93000 

1.31800 

1.39200 

1.46600 

1.54000 

2.38000 

1.56000 

.74000 

-.08000 

^.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.08034 

.10276 

.13592 

.17514 

2.18000 

2.17000 

2.16000 

2.15000 

1.39500 

1.45700 

1.51800 

1.58000 

.83000 

.54000 

.25000 

-.04000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- .06933 

— .08792 

.09070 

.11225 

2.39000 

2.38700 

2. 38300 

2.38000 

- 1.44500 

1.49700 

1.54800 

1.60000 

. 13000 

.09700 

.06300 

.03000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.04704 

.05687 

.06238 

.06696 

2.59000 

2.58700 

2.58300 

2.58000 

1.46800 

1.51200 

1.55600 

1.60000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.03515 

.03535 

.03554 

.03601 

2.77430 

2.77430 

2.77430 

2.77430 

1.47890 

1.51930 

1.55960 

1 .60000 

0.00000 

0.00000 

0.00000 

0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.03010 

.03008 

.03008 

.03009 

-- 2.92390 

- 2.92390 

2.92390 

2.92390 

1.47890 

1.51930 

1.55960 

1.60000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

.03596 

.03596 

.03596 

.03596 
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TABLE II . ROTOR FLOW PATH AND BLADE GEOMETRY 
FOR /3g = 10.0 DEGREES 


Rotational Speed = 95,000 rpm 
Inducer Tip Speed = 2250 fps 

STATION NUMBfS 1 


, p 

1.7b300 

1 .61A17 

1.A7533 

1.33650 

1.19767 

1.05883 

.92000 

1 

-.20000 

-.20000 

-.20000 

-.20000 

-.20000 

-.20000 

-.20000 

^ BETA 

60.50000 

57.70000 

5A. 90000 

52.10000 

A9. 01700 

A5. 93300 

A2. 85000 

,0EV 

-0.00000 

-0.00000 

-0.00000 

-O.OCCOO 

-0 .00000 

-0.00000 

-0.00000 

THICK 

.02781 

.02781 

.03532 

.0A771 

. 06AA8 

.083A1 

. 10701 

STATION 

NUMBER 2 
1.75300 1.61A17 

1 .A7533 

1. 33650 

1 .19767 

1 .05883 

.920C0 


-.13330 

-. 13330 

-.13330 

-.13330 

-.13330 

-.13300 

-.13300 


6Q.A5000 

57 .53300 

5 A. 61700 

5 1. 7C0C0 

A8 .6 6300 

A5 .66700 

A? .65000 

• DEV 

-0.00000 

-0. 00000 

-0 .00000 

-0. OOOCO 

-0.00000 

-0.00000 

-0. OOOCO 

.THICK 

.02781 

. 02 781 

.03532 

.0A771 

.06AAB 

. C83A1 

.10701 

-—STATION 

‘P 

NUMBER 3 
1.75300 1.61A80 

1 .A7670 

1.33850 

1. 2C030 

1 .06220 

.92AOO 

-,2 

-.06670 

-.06670 

-. C6670 

-.06670 

-.06670 

-.C6670 

-.06670 

BETA 

60. 10000 

57.13300 

5 A. 16 700 

51.20000 

A8. 20000 

A5 .20000 

A2. 20000 

^DEV 

-0.00000 

-0. 00000 

-0 .00000 

-0.00000 

-O.OCCOO 

-0.00000 

-0, CCCCO 

"THICK 

1 

.02781 

.02781 

. 03532 

.0A771 

. 06AA8 

.C83A1 

.10701 

STATION 

NUMBER A 
1,75300 1.61530 

1.A7770 

1 .3A000 

1. 2C230 

1 .C6A70 

.92700 

i Z 

0. OCOOO 

0 .00000 

0 . oooco 

O.OCCOO 

0 .00000 

0, CCCCO 

0 .00000 

. beta 

59. AOOOO 

56.A3300 

53 .A6700 

50.50000 

A7. 50000 

AA .50000 

Al. 5C000 

DEV 

-0.00000 

-0.00000 

-C. 00000 

-0.00000 

-0. OCCCO 

-O.CCOCO 

-0.00000 

^THICK 

. 02781 

.02781 

.03532 

.0A771 

.06AA8 

.C83A1 

.09865 

STATION 

R 

NUMBER 5 
1.75300 1.61720 

1.A8130 

1.3A550 

1.2C970 

1 .07360 

.93800 

■”Z 

.06700 

.07550 

.08A00 

.09250 

. 10100 

.10950 

.11800 

luBETA 

58.A5000 

55.30000 

52.15000 

A9. 00000 

A6 .00000 

A3. 00000 

AO. 00000 

jiDEV 

-0.00000 

-0.00000 

-0,00000 

-0.00000 

-O.OCCOO 

-0.00000 

-0.00000 

,, THICK 

.028A2 

.028A2 

.03532 

.0A82A 

.06513 

,08A75 

.1075A 

STATION 

mP 

NUMBER 6 
1.75300 1.61920 

1. A8530 

1.35150 

1.21770 

1.08380 

.95000 

4 SrT 

.13600 

.15A50 

. 17300 

.19150 

.21000 

.22850 

.2A700 

beta. 

56..90000 53.66700 - 

50.A3300 

A7. 20000 - 

A3. 96700 

AO. 73300 

37. 5C000 

1. DEV 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-O.OCCOO 

-C.COOOO 

-0.00000 

4, {thick 

.02861 

.02901 

.03573 

.0A830 

.06592 

.08691 

.10800 

■^1 _ station 

NUMBER 7 
1.75300 1.62250 

1.A9200 

1.36150 

1.23100 

1.10050 

.97000 

z 

.21000 

.23800 

,26700 

.29500 

.32300 

.35170 

.38000 

BETA 

5A. 67000 

51.28000 

A7, 89000 

AA.A9000 

Al.lCCCO 

37.71000 

3A. 32000 

-'DEV 

-0.00000 

-0.00000 

-0.00000 

.03586 

-0.00000 

.0A878 

-0.00000 

.06707 

-C. 00000 
.08810 

-0.00000 

.11517 

. thick 

1 

.02875 

.02965 

"'1 station 

'-'‘P _ 

NUMBER 8 
1.75300 1.62750 

1.50200 

1.37650 

1.25100 

1.12550 

1.00000 


.28500 

.32A50 

.36A00 

. A0350 

. AA300 

. A8250 

. 52200 

„BETA 

52.00000 

A8. AlOOO 

AA. 83000 

A1 .2A000 

37.65C00 

3A. 07000 

30. AROOO 

DEV 

'c'lHlCK 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0,00000 

-0.00000 

-0.00000 

.D2B85_ 

__.03002 „ 

. .036A5 

- .0A937- 

.0673A 

.08860 

.12580 

STATION 

R 

NUMBER 9 
1.75A00 1.63670 

1.51900 

1. AO2C0 

1 .28500 

1. 16700 

1.05000 

- 2 

.36000 

. A1500 

. A6900 

.52A00 

. 57900 

.63300 

.68800 

BETA 

A8 .8 5000 

AA. 96000 

A1.070CC 

37.17000 

33.28000 

29 . 390C0 

25.50000 

QEV 

-0.00000 '-0. 00000"' 

-0.00000 

-0.00000 " 

-0.00000 

■ -0.00000 

-0. CCOOO 

|THIC K 

.02901 

.03057 

.03679 

.0A9A3 

. 06772 

.09A75 

. 1265A 

' 1 STATION 

P, 

NUMBER 10 
1.76000 1.65A00 

1 .5A800 

1. AA3C0 

1.33700 

1 .23100 

1. 125C0 

|Z 

■BETA 
fTo EV 

.AA500 

. 5 1500 

. 58 ACC 

.65A00 

. 72300 

.79300 

. 86500 

AA. 78000 . AO. 60000 
-0.00000 -0.00000 

36.A3000 
-0 .00000 

32.25000 
-0. OOOCO 

28.07000 

-0.00000 

23.9CCC0 

-0.00000 

19.72000 
-0. OOOCO 

Thick 

.02927 

.03080 

.03798 

.05086 

.07156 

.09593 

.1A22A 



Station Number 4 = 

Blade Trailing 

Edge 


Station Number 19 = 

Blade 

Leading 

Edge 

STATION 

NUMBER 11 

1 





‘ R 

1.77500 1.6BAC0 

1.59300 

1.50300 

1. A12C0 

1 .32100 

1. 230C0 

2 

.5A500 .62700 

. 709 CO 

.79100 

.87300 

.9 5500 

1 .OAOOO 

BETA 

39.3AOOO 35.09000 

30,85000 

26.6C0CC 

?2 .35000 

18.11000 

13.86000 

' DEV 

-0.00000 -0.00000 

-0 .00000 

-0.00000 

-0. CCCOO 

-0 .00000 

-0.00000 

' THICK 

.02931 .03107 

1 .03803 

.05136 

.07159 

.09816 

.151A8 

’ ' STATION 

NUMBER 12 






' R 

1.80000 1.72600 

1.65300 

1.57900 

1. 5C500 

1.A3200 

1.35800 

M 2 

.65500 .7A700 

.8 38 CC 

.93000 

1 .02200 

1.11300 

1 .20500 

,, B ETA 

32.36000 28.35000 

2A .33000 

20.320CC 

16.31000 

12.29000 

8 . 28000 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0. OCCOO 

-0 .OCOOO 

-0.00000 

THICK 

- .029A7 .0312A 

. 03816 

.05302 

.07266 

. 1C703 

.15825 

',4! STATION 

NUMBER 13 






R 

1.8A000 1.78600 

1.73100 

1.67700 

1.62200 

1.56800 

1.51300 

'■ 2 — 

— .78700 - .88300 

-- .97800 

1.07A00 

1 .16900 

1.265C0 

1 .35780 

1 . BETA 

23.A8000 20.03000 

16.58000 

13.130CC 

9 .68C00 

6 .23000 

2. 780C0 

I'DEV^^ 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-C. OCOOO 

-0.00000 

THICK 

II 

.02952 .03180 

.OAO A6 

.05711 

.0E62 7 

.12351 

.17630 

m! station 

'number 1A 






41. E 

1.92500 1.88750 

1.850C0 

1.81250 

1. 77500 

1.73750 

1 .70000 

2 

.9AOOO 1,02800 

1. 11500 

1.20300 

1 .29000 

1.37800 

1. A6500 

'■BETA- 

-13. 59000--11. 03000 

--*.a.A7000 

- 5,91000 

3.35000 

.79000 

-1.77000 

-DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0. OCOOO 

-0.00000 

-0.00000 

JHICK 

,02981 .03A5A 

.0A758 

.06803 

.10223 

.13A18 

.183AS 

'M - STATION 

NUMBER 15 






4 R 

2.05600 2.03500 

2.01A00 

1 .99300 

1.97200 

1.95100 

1.93000 

>1 2 

1.09500 1.16900 

1.2A300 

1.31800 

1.39200 

1.A6600 

1 .5A000 

beta 

A. 20000 2,50000 

.80000 

-.90000 

-2.6C000 

-A. 30000 

-6,00000 

-'DEV 

_5-0, 00000 --0. 00000 

^0.00000 

-0.00000 

-0.00000 

-0. COOOO 

-0.00000 

ijTHlCK 

.03299 .0A162 

.05690 

.0803A 

.10276 

.13592 

.1751A 

■'I STATION 

NUMBER 16 






“R -- 

2.21000 2.20000 

2.19000 

2,18000 

2.17000 

2.16000 

2.15000 

'2 ^ 

1.21000 1.27200 

1,33300 

1,39500 

1.A5700 

1.51800 

1.58000 

1 BETA 

-2.50000 -3.A8000 

-A, A5000 

-5.A3000 

-6. ACOOO 

-7.38000 

-B. 35000 

, DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

' THICK— 

VI I 

- -.03515— .0A363 

.05579 

- .06933 

.08792 

.09870 

.11225 

I I STATION 

NUMBER 17 







2.A0000 2.39700 

2.39300 

2.39000 

2 .38700 

2.38300 

2.38000 

2 

- 1.29000 1.3A200 

1.39300 

1.AA500 

1 . A9700 

1.5A800 

1.60000 

4'BETA 

-7.20000 -7.63000 

-8.05000 

-8.A8000 

-8.90000 

-9.33000 

-9.75000 

4 ., DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

, THICK 

.03A82 .0379A 

1 .0A252 

.0A70A 

.05687 

.06238 

.06696 

■i ' STATfON 

“number 13 






„R 

2.60000 2.59700 

2.59300 

2.59000 

2 .58700 

2.58300 

2.58000 

I L, 

1.33500 1.37900 

1 .A2300 

1.A6800 

1.51200 

1.55600 

1.60000 

‘ BETA 

-9.70000 -9.75000 

-9.80000 

-9.85000 

-9.90000 

-9.95000 

-10.00000 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0 .00000 

-0.00000 

-0. COOOO 

1 THICK 

1 

.03A71 .03A86 

.03A98 

,03515 

. 03535 

.0355A 

.03601 

— STATION-NUHBER 19- 






"R 

2.71A00 2.71A00 

2.71AOO 

2.71A00 

2.71A00 

'2.71A00 

2.71A00 

•*/ z _ 

1.35120 1.39270 

1. A3A10 

— 1.A7560 

1.51710 

1.55850 

1 .60000 

BETA -10.00000 -10.00000 

-10.00000 

-10.00000 

-10.00000 

-10.00000 

-10.00000 

DEV 

-0.00000 -0.00000 

-0.00000 

-0,00000 

-0. CCCOO 

-0.00000 

-0. OOOCO 

'THICK 

.03038 .03025 

I . 03013 

.03010 

.03008 

.C3008 

.03009 

STATION 

NOTW^ 20 






-R 

— 2.86780 2, 86780 

- 2.86780 

2.86780 

2.86780 

2.8 67 80 

2 . 86780 

2 

1.35120- 1 .39270 

1. A3A10 

1.A7560 

1 .51710 

1.55850 

1 .60000 

1 , BETA 

-10.00000 -10.00000 

-10.00000 

-10. OOOCO 

-10.00000 

-10.00000 

-10. COOOO 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0. OCCOO 

-C. CCOOO 

-0.00000 

'THICK 

- .03522 .03522 

. 03522 

.C3522 

.03522 

.02522 

.03522 


137 



This Page Intentionally Left Blank 



TABLE III. ROTOR FLOW PATH AND BLADE GEOMETRY 
FOR /3g = 20.0 DEGREES 


Rotational Speed = 95, 

000 rpm 







Station Number 5 

= Blade 

Trailing 

Edge 

Inducer Tip Speed = 2200 fps 







Station Number 20 

= Blade 

j Leading 

Edge 








STATION 

NUMBER 11 


1 












R 

1,76000 

1.65A00 

1.5A800 

1.AA300 

1.33700 

1.23100 

1.12500 

STATiaN 

NUHPER 1 






.Z 

. AA500 

.51500 

. 58A00 

.65A00 

.72300 

. 79300 

.86500 

R 

1.75300 1.61A17 

1.A7533 

1.33650 

1.19767 

1.05883 

.92000 

BETA 

39.50000 

35,22000 

30.93000 

26.65000 

22.37000 

18.08000 

13.30000 

z 

-.30000 -.30000 

-.30000 

-.30000 

-.30000 

-.30000 

-.30000 

DEV 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0. 00000 

-0.00000 

BETA 

60.60000 58.5<tOOO 

56.70000 

53.53000 

50. A8600 

A6. 99000 

A2. 96600 

THICK 

.02927 

. 03080 

.03798 

.05086 

.07156 

.09593 

.1A22A 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 









THICK 

.02781 .02781 

.03532 

.0A771 

.06AA8 

.083A1 

.10701 

STATION 

NUMBER 12 














R 

1.77500 

1.68A00 

1.59300 

1,50300 

1. A1200 

1. 32100 

1.23000 

STATION 

NUMBER 2 





1 

Z 

. 5A500 

.62700 

1 .70900 

.79100 

.87300 

.95500 

l.OAOOO 

R 

1.75300 1.61A17 

1.A7533 

1.33650 

1.19767 

1.05883 

.92000 

BET A 

33.00000 

28.51000 

2A. 02000 

19,52000 

15.03000 

10. 5AOOO 

6.05000 

Z 

-.20000 -.20000 

-.20000 

-.20000 

-.20000 

-.20000 

-.20000 1 

DEV 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

60.50000 57,70000 

5 A. 90000 

52.10000 

A9. 01700 

A5. 93300 

A2. 85000 

THICK 

.02931 

.03107 

.03803 

.05136 

.07159 

.09816 

.151 A8 

DEV 

-O.COCCO -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 









THICK 

.02781 .02781 

.03532 

.OA771 

.06AAB 

.083A1 

.10701 

STATION 

NUMBER 13 














R 

1 .80000 

1.72600 

1.65300 

1.57900 

1.50500 

1. A3200 

1.35800 

STATION 

NUMBER 3 





i 

Z 

. 65500 

.7A700 

.83800 

.93000 

1.02200 

1.11300 

1.20500 

R 

1.75300 1.61'.17 

1 .A7533 

1.33650 

1.19767 

1.05883 

.92000 

BET A 

25.75000 

21.33000 

16.90000 

12.A8000 

8.05000 

3. 63000 

-.80000 

Z 

-.13330 -.13330 

-.13330 

-.13330 

-.13330 

-.13300 

-.13300 

DEV 

-0, 00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

50.80000 56.88000 

53.97000 

51.05000 

A8. 13000 

A5. 22000 

A2. 30000 

THICK 

.029A7 

.0312A 

.03816 

.05302 

.07266 

. 10703 

.15825 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 









THICK 

.02781 .02781 

.03532 

.0A771 

.06AA8 

.083A1 

.10701 

STATION 

NUMBER lA 














R 

1.8A000 

1.78600 

1.73100 

1.67700 

1.62200 

1.56800 

1.51300 

STATION 

NUMBER A 






z 

. 78700 

.88000 

.97310 

1.06610 

1.15910 

1.25220 

1.3A520 

R 

1.75300 1.61A80 

1.A7670 

1.33850 

1.20030 

1.06220 

.92A00 1 

BETA ~ 

16.20000 

12.A6000 

8.73000 

A. 99000 

1.25000 

-2. A8000 

-6.22000 

Z 

-.06670 -.06670 

-.06670 

-.06670 

-.06670 

-.06670 

-.06670 

OEV 

-0.00000 

-0.00000 

-0.00000 

-0,00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

59.18000 56.77000 

53.35000 

50. AAOOO 

A7. 53000 

AA. 61000 

Al. 70000 

THICK 

. 02952 

.03180 

.0A0A6 

.05711 

.08627 

.12351 

.17630 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 ' 




1 





THICK 

.02781 .02781 

.03532 

.0A771 

.06AA8 

.083A1 

.10701 , 

STATION 

NUMBER 15 














R 

1,92500 

1.88750 

1.85000 

1.81250 

1.77500 

1.73750 

1.70000 

STATION 

NUMBER 5 






Z 

.9AC00 

1.02800 

11.11500 

1.20300 

1.29000 

1.37800 

1.A6500 

R 

1.75300 1.61530 

1. A7770 

1.3A000 

1.20230 

1.06A70 

.92700 

BETA 

5. AAOOO 

2.62000 

-.19000 

-3.01000 

-5.83000 

-8.6A000 -11.A6000 

Z -- - - 

0.00000 0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 1 

OEV — 

-0.00000 

-0.00000 

- -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

58.20000 55.30000 

52.AQ000 

A9. 50000 

A6. 63000 

A3. 76700 

AO. 90000 ! 

THICK 

.02981 

.03A5A 

1 .0A758 

.06803 

.10223 

.13A18 

.183A5 

DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 1 









THICK 

.02781 .02781 

.03532 

.0A771 

.06AA8 

.083A1 

.10701 ' 

STATION 

NUMBER 16 













1 

R 

2.05600 

2.03500 

2.01A00 

1.99300 

1.97200 

1.95100 

1.93000 

STATION 

NUMBER 6 






Z 

1.09500 

1.16900 

1.2A300 

1.31800 

1.39200 

1. A6600 

1.5A000 

R 

1.75300 1.61720 

1.A8130 

1. 3A550 

1.20970 

1.07380 

.93800 

BET A 

-6.00000 

-7.60000 

-9.20000 

-10.80000 

-12.AOOOO 

-lA. 00000 - 

15.60000 

Z 

.06700 .07550 

.08A00 

.09250 

.10100 

.10950 

.11800 i 

OEV 

-0.00000 

-0. 00000 

-0. 00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

1 BETA 

55.50000 57.75000 

50.00000 

A7. 25000 

AA. 50000 

Al. 75000 

39.00000 j 

THICK 

- .03299 - 

- .0A162 

.05690 

.0803A 

.10276 

.13592 

.1751A 

1 DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0. 00000 

-0.00000 ' 

, 








1 THICK 

.028A2 .028A2 

.03532 

.0A82A 

.06513 

.08A75 

.1075A 1 

STATION 

NUMBER 17 













1 

R 

2.21000 

2.20000 

2.19000 

2.18000 

2.17000 

2.16000 

2.15000 

STATION 

NUMBER 7 





1 

Z 

1.21000 

1.27200 

1.33300 

1.39500 

1. A5700 

1.51800 

1.58000 

R 

1.75300 1.61920 

1.A8530 

1.35150 

1.21770 

1.08380 

.95000 1 

BETA 

-13.30000 - 

-lA. 12000 

-lA. 93000 

-15.75000 

-16.57000 

-17.38000 - 

■18.20000 

Z 

.13600 .15A50 

.17300 

. 19150 

.21000 

.22850 

.2A700 i 

OEV 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

53.50000 50.55000 

A7. 60000 

AA. 65000 

Al. 70000 

38.75000 

35.80000 ■ 

THICK 

.03515 

.0A363 

.05579 

.06933 

.08792 

.09870 

.11225 

1 DEV - - 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 i 

— 


- 






THICK 

.02861 .02901 

.03573 

.0AB30 

.06592 

.08691 

.10800 

STATION 

NUMBER 18 














' i R 

2. AOCOO 

2.39700 

2.39300 

2.39000 

2.38700 

2. 38300 

2.38000 

STATION 

NUMBER 8 







1.29000 

1.3A200 

1.39300 

I.AA500 

1. A9700 

1.5A800 

1.60000 

R 

1.75300 1.62750 

1. A9200 

1.36150 

1.23100 

1.10050 

.97000 ' 

BETA -18.00000 -18.30000 

-18.60000 

-18.90000 

-19.20000 

-19.50000 -19.80000 

Z 

.21000 .23800 

.26700 

.29500 

.32300 

.35170 

.38000 

DEV 

-0.00000 

-0.00000 

-0.00000 

-0,00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

50.55000 A7.A6000 

AA. 37000 

Al. 27000 

38.18000 

35.09000 

32.00000 

THICK 

. 03AB2 

.0379A 

.0A252 

.OA70A 

.05687 

.06238 

.06696 

, DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 1 









^HICK 

r02875 -.02965- 

.03586 

.0A878 

-- .06707 

— .08810 

.ii5ir 1 

, -STATION' 

-NUMBER 19 


1 











1 

' R 

2.60000 

2.59700 

2.59300 

2.59000 

2.58700 

2.58300 

2,58000 

STATION 

NUMBER 9 





1 

! Z 

1. 33500 

1,37900 

1. A2300 

1.A6800 

1.51200 

1.55600 

1.60000 

R 

1.75300 1.62750 

1.50200 

1.37650 

1.25100 

1.12550 

1.00000 

BtIA -20.00000 -20.00000 

-20.00000 

-20.00000 

-20,00000 

-20.00000 -20.00000 

~ ■ 

Z 

.28500 .32A50 

.36A00 

.A0350 

.AA300 

.A8250 

.52200 

UfcV 

-O.OUOOO 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

BETA 

-A7»70000 - AA. 78000 

-A0.87000— 37.A5000 

- 3A. 03000 

30.62000 - 

27.20000 

THICK 

.03A71 

.03A86 

.03A98 

.03515 

.03535 

.0355A 

.03601 


DEV 

-0.00000 -0.00000 

-0.00000 

-0,00000 

-0.00000 

-0.00000 

-0.00000 










THICK 

.02885 .03002 

.036A5 

.0A937 

.0673A 

,08860 

.12580 

STATION 

NUMBER 20 


1 












1 

ff — 

- 2.65370 

- 2,65370 

2.65370 

2.65370 

2.65370 

2.65370 

2.65370 


STATION 

NUMBER 10 





1 

i 1 \ 

1. 3A520 

1.38770 

11. A3010 

1.A7260 

1.51510 

1.55750 

1.60000 


R 

1.75A00 1.63670 

1.51900 

1.A0200 

1.28500 

1. 16700 

1.05000 ! 

BtIA -ZO. 00000 -20.00000 

-20.00000 

-20.00000 

-20.00000 

-20.00000 -20.00000 


1 Z 

.36000 .A1500 

.A6900 

.52A00 

,57900 

.63300 

.68800 

utv 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0. 00000 

-0.00000 


BETA - - 

A5.3A000 Al. 25000 

37.16000 

33.07000 

28.98000 

2A. 89000 

20.80000 

THI CK 

.03038 

.03025 

1 .03013 

.03010 

.03008 

.03008 

.03009 

1 DEV 

-0.00000 -0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 









THICK 

.02901 .03057 

.03679 

.0A9A3 

.06772 

.09A75 

.1265A 

STATION 

NUMBER 21 













1 

1 5 

2.81110 

2.81110 

2.81110 

2.81110 

2.81110 

2.81110 

2.81110 








V 2 

- 1.3A520 

- 1.38770 

1. A3010 

1.A7260 

1.51510 

1.55750 

1.60000 








BETA 

-20.COOOO -20.00000 

-20.00000 

-20.00000 

-20.00000 

-20.00000 - 

20.00000 








DEV 

-0. 00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 








THICK 

. 03A58 

.03A58 

, .03A58 

.03A58 

.03A58 

.03A58 

.03A58 
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Rotational Speed = 95,000 rpm 
Inducer Tip Speed = 2250 fps 

NOTE; 


TABLE IV. OPTIMIZED TANGENTIAL THICKNESSES 
FOR /3g = 10 DEGREES 


X = Axial distance (inches) 
T = Tangential thickness (i 

S TIT I DM NO. 1 


STREAMLINE X 

R 

T 

1 

• 2468 

2.7140 

.0417 

2 

. 2261 

2.7140 

.0427 

3 

.2073 

2.7140 

.0433 

k 

• 1666 

2.7140 

.0492 

5 

• 1669 

2.7139 

.0494 

6 

• 1451 

2.7140 

• 0496 

7 

• 1244 

2.7140 

.0548 

8 

• 1037 

2.7140 

.0539 

9 

• 0629 

2.7140 

.0542 

10 

• 0622 

2.7140 

.0575 

11 

• 0415 

2.7140 

.0560 

12 

.0207 

2.7140 

• 0560 

13 

0.0000 

2.7140 

.0587 


STATION 

NO. 2 



STREAMLINE X 

R 

T 

1 

.2613 

2.6235 

• 0484 

2 

• 2388 

2.6226 

• 0541 

3 

• 2163 

2.6223 

• 0562 

4 

• 1938 

2.6217 

• 0568 

5 

.1713 

2.6211 

.0637 

6 

• 1469 

2.6205 

.0659 

7 

.1264 

2.6199 

.0699 

B 

• 1039 

2.6193 

.0705 

9 5 

• 0814 

2*6167 

.0718 

10 

• 0569 

2.6181 

.0746 

11 

.0364 

2.6175 

.0736 

12 

• 0139 

2.6169 

.0743 

13 

-.0086 

2.6163 

.0748 


station 

NO. i 



streamline X 

R 

T 

1 

.2771 

2.5331 

• 0565 

2 

• 2525 

2.5318 

• 0628 

3 

• 2280 

2.5306 

• 06 60 

4 

• 2035 

2.5294 

.0704 

5 

. 1790 

2.5282 

.0776 

6 

.1544 

2.5270 

.0822 

7 

.1299 

2.5258 

.0857 

8 

.1054 

2.5246 

.0908 

9 

• 0609 

2.5234 

.0917 

10 

.0563 

2.5221 

.0925 

11 

• 0318 

2.5210 

.0950 

12 

.0073 

2.5198 

.0956 

13 

-.0172 

2.5185 

.0958 


STATION NO. A 

STREANLINE * P T 


1 

.2977 

2.4435 

.0621 

8 

• 1326 

2.1632 

• 1523 






nches) z 

.2707 

2.4415 

.0668 

9 

.0939 

2.14B3 

.1669 


station 

NO. 11 



3 

. 2437 

2.4397 

.0746 

10 

.0552 

2.1434 

.1746 






4 

• 2166 

2.4378 

• 0823 

11 

.0165 

2.1386 

.1779 


streamline X 

R 

T 

5 

. 1896 

2.4359 

• 0894 

12 

-.0222 

2.1336 

.1835 






6 

• 1626 

2.4340 

.0979 

13 

-.0609 

2.1287 

.1873 

", 

1 

.6600 

1.9249 

.0535 

1 1 

• 1356 

2.4322 

.1007 





’ 

2 

.5693 

1.9084 

• 0614 

' 8 

. 1086 

2.4303 

.1092 





I i 

3 

• 5185 

1.8918 

• 0693 

L] 9 

• 0816 

2.4284 

• 1106 

STATION 

NO. 8 




4 

.4478 

1.8752 

• 0832 

10 

• 0546 

2.4265 

.1145 





! 1 

5 

.3770 

1.8566 

• 0985 

11 

.0276 

2.4247 

• 1148 

STREANLINE X 

R 

T 

1 

6 

• 3063 

1.3421 

• 1264 

12 

• 0006 

2.4228 

• 1134 






7 

• 2355 

1.8254 

• 1803 

13 

-•0264 

2.4209 

• 1151 

1 

.4572 

2.1116 

.0607 


8 

• 1648 

1.8009 

• 2180 





2 

.4127 

2.1049 

.0679 


9 

.0941 

1.7923 

• 2441 





3 

• 3682 

2.0983 

.0776 


10 

.0233 

1.7757 

.2743 

M STATION 

NO. 5 



4 

. 3237 

2*0916 

.0869 

! 

11 

-.0474 

1.7591 

• 3051 





5 

.2792 

2.0850 

• 1084 

' 1 

12 

-.1046 

1.7594 

.3282 

STREAMLINE X 

R 

T 

6 

.2346 

2.0703 

.1348 


13 

-.1600 

1.7500 

• 3424 





7 

.1903 

2.0717 

• 1436 






' 1 

• 3248 

2.3552 

.0627 

8 

• 1458 

2.0650 

.1672 






' 2 

.2947 

■'2.3525 

.0741 

9 

• 1013 

2.0594 

• 1850 


STATION 

NO. 12 “ 



3 

• 2646 

2.3498 

.0775 

10 

• 0568 

2.0516 

• 1946 






4 

• 2345 

2.3472 

• 0936 

11 

• 0123 

2*0451 

• 2046 


STREAMLINE X 

R 

T 

5 

.2044 

2.3445 

.0976 

12 

-.0322 

2.0385 

.2146 






u 6 

.1743 

2.3419 

.1137 

13 

-.0767 

2.0316 

• 2223 

‘ , 

1 

.7730 

1.8584 

• 0519 

. . 7 

• 1442 

2.3393 

.1171 





i , 

2 

.7033 

1.8390 

.0573 

' '■ 8 

• 1141 

' 2.3366 

.1279 





1 ' 

3 

• 6337 

1.8195 

• 0616 

,1 9 

• 0841 

2.3340 

.1299 

STATION 

NO. 9 



1 

4 

.5641 

1.8001 

• 0716 

' 10 

• 0540 

2.3313 

• 1344 





t 

5 

• 4945 

1.7808 

.0805 

11 

• 0239 

2.3287 

.1350 

streamline X 

R 

T 


6 

.4248 

1.7613 

• 1011 

12 

-.0062 

2.3260 

.1355 






7 

• 3552 

1.7419 

.1290 

13 

-•0363 

2.3233 

.1357 

1 

.5187 

2.0419 

.0591 


8 

• 2856 

1.7225 

.1733 

1 




2 

.4673 

2.0331 

.0649 


9 

.2159 

1.7031 

.2149 

' 




3 

• 4160 

2.0243 

.0734 


10 

• 1463 

1.6637 

.2500 

, STATION 

NO. 6 



4 

.3647 

2.0154 

• 0865 

, 

11 

.0767 

1.6643 

• 2822 





5 

• 3133 

2.0066 

• 1045 






STREAMLINE X 
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Rotational Speed = 95,000 rpm 
Inducer Tip Speed = 2200 fps 
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